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ABSTRACT 

Context. Abridged 

Aims. The goal of this paper is to quantify the changes on the SF distribution within the disc galaxies in the last ~8 Gyr We use as a 
proxy for the SF radial profile the Near-UV surface brightness distributions, allowing suitably for extinction. 

Methods. We compare the effective radii (Reff ) and concentration of the flux distribution in the rest-frame Near-UV for a sample of 
270 galaxies in the range < z < 1 . This radial distribution is compared to that measured in the rest-frame B -band, which traces older 
stellar populations. The analysis is performed using deep, high resolution, multi-band images from GALEX , SDSS , and HST/ACS 

- GOODS -South. 

Results. The relation R^ff (NUV)-M-i, suffers a moderate change between z ~1 and z ~0, at a fixed stellar mass of 10'" Mq galaxies 
increase their effective radii by a factor 1.18±0.06. The ratio R^ff (NUV )/Reff (S) has increased by ~10% in the same period. Median 
profiles in A'6' V show signs of truncation at R ~Reff , and median colour profiles (NUV -B) show a minimum (a "bluest" point) also 
around R~l-l.5Reff ■ The surface brightness of the discs (at R^ff ) has decreased by ~80% in NUV , and by ~60% in 6 since z ~1. 
The distributions of NUV flux are more compact at z~l than nowadays, in terms of the fraction of flux enclosed in a specific radius 
(in kpc). The central brightness in i?and NUVhas increased with respect to the brightness level of the discs, both in Sand NUV , 
indicating perhaps a continuous accumulation of the bulge structure since z ~1, or that they have been left exposed because of the 
significant decline in the the star-forming activity in the discs. 

Conclusions. Our results indicate that the SF surface density has decreased dramatically in discs since z ~1, and this decline has been 
more intense in the central parts (<Reff ) of the galaxies. In addition, our data suggest that the bulges/pseudo-bulges have grown in 
surface brightness with regard to the discs since z~l. 

Key words. Galaxies: evolution - Galaxies- stellar content - Galaxies: high-redshfit - Galaxies: photometry 



' 1. Introduction 

Observations in the last decade have shown that the global star 
formation rate (SFR) in galaxies has suffered a significant de- 
' crease, by roughly tenfold, since z~l (e.g. Lilly et al. 119961 
Madau et al. 119961 11998I I. Other observations have shown that 
this star-forming activity has gradually shifted to less massive 
galaxies, a process commonly termed as "downsizing" (e.g. 
Cowie et al. 119961 Brinchmann & Ellis 2000). This dramatic 
change in the pace at which star formation (SF hereafter) has 
taken place, and the shift in the masses of the hosts which har- 
bour this activity are yet to be fully understood. Also impor- 
tant, and pending, is to find a consistent scheme in which this 
evidence is related to the evolution in the mass content (e.g. 
Rudnick et al. |2003l Bell et al. |2007l Perez-Gonzalez et al . 
12008 ), lumino sities ( e.g. Ilbert et al. 2006 Conselice et al. 120051 
Rudnick et al. 120031 Wolf et al. 2003.) a nd stru cture (morpholo- 
gies) of galaxies (e.g. Conselice et al. 120081 Coe et al. 120061 
Cassata et al.'200?, Conselice |2003l l. 

An area of study which has not yet received enough atten- 
tion, and could provide clarifying evidence to solve the prob- 
lems posed, is that of the spatial distribution of the SF within 



galaxies as a function of redshift. In particular, we refer to stud- 
ies of the radial distribution of SF as a function of epoch. This 
kind of work adds constraints on the process of galaxy evolu- 
tion, as knowing where stars are being born in galaxies during 
the history of the universe is essential to reconstruct their evo- 
lutionary history. This article is directed at this research field, 
surveying the radial distribution of Near-UV [NUV ) flux, used 
as a tracer of SF, in a sample of disc galaxies within the redshift 
range 0<z<L 

There are several pieces of evidence which suggest that 
undisturbed, i.e. non interacting, disc galaxies are the most 
promising subjects to study the change of the spatial distribution 
of SF with cosmic time. First of all, these are the most numerous 
galaxies in the Local Universe, accounting for ~70% of the to- 
tal, and their preponderance in number density seems to hold, at 
least, up to z ~ 1 (Conselice et al. 2005 ). Moreover these galaxies 
harbour most of the unobscured SF at z <1 (Bell et al. 120051) . It 
has been hypothesized that the decline in global SFR between 
z ~l and z ~0 is due to a decreased rate of major mergers (e.g. 
Somerville et al. 2001), but there is evidence to challenge this 
picture. For example. Bell et al. ( 120051 ) found that more than half 
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of the intensely star-forming galaxies at z~0.7 have spiral mor- 
phologies, without major disturbances, whereas less than 30% 
are strongly interacting. They proposed alternative explanations 
for the decrease in global SFR, involving gas depletion and weak 
interactions with small satellites. In a complementary line of re- 
search. Bell et al. 120071 explored the relation between star for- 
mation and the acquisition of stellar mass for galaxies at z <1. 
They found that the increase in stellar mass content in this time 
interval is consistent with the integrated SFR. Most of these stars 
are formed in galaxies located in the "blue cloud" (typically late- 
type galaxies with young stellar populations), but there must be 
some mechanism of star formation quenching by which part of 
these galaxies pass to the "red sequence" (galaxies which are 
passively evolving, with little or no star-forming activity, though 
there are cases in which they harbour significant dust-obscured 
SF, and which usually have early morphological types) at some 
point in that period. Otherwise, there would be a dramatic over- 
production of blue cloud galaxies in the local universe over the 
observed populations. For these reasons, it is clear that undis- 
turbed disc galaxies are an interesting target for a study of the 
spatial distribution of SF. 

An aspect of spiral galaxies which has received much atten- 
tion in recent years is the formation and evolution of their stellar 
discs. Recent studies of galaxies at z >1 in the HS T -Ultra Deep 
Field have shown that most late-type galaxies at those times 
present morphologies which are dominated by massive clumps 
of stars (10^ <M* < lO** Mq ), which some authors argue can dis- 
solve, by dynamical friction, into exponential discs in time scales 
of the order of 1 Gyr (see e.g. Bournaud et al. 2007; Bournaud 
et al. 2008; Elmegreen et al. 2008a,b). Between z =1 and z -Q> 
the stellar discs, which show more regular morphologies, have 
grown in effective radius, R^ff , measured in B -band, so as to 
keep their surface mass densities roughly constant (Barden et 
al. [2005). This means, that there is no evolution in the Reff- 
Mi, relation in that period, although each individual galaxy must 
evolve along this relation. In addition, from the study of stellar 
disc truncations in this range of redshifts, it seems that there is a 
moderate growth, by a 25-30%, in the radii of these truncations 
at a given stellar mass of the galaxies (Trujillo & Pohlen 2005] 
AzzoUini et al. l2008bl l. Moreover, there is growing evidence that 
supports the idea that this growth, at least in relatively recent 
times, has been inside-out. Observing how this apparent growth 
of the stellar discs actually takes place, by building a history of 
the spatial distribution of SF with epoch is vital to reach a sound 
knowledge of the underlying processes. 

The study of the spatial distribution of SF has been so far 
restricted to galaxies in the nearby universe. The early work of 
Paul Hodge and collaborators was pioneering in this field (e.g. 
Hodge [T#9al[T969b; Hodge & Kennicutt 1983 ). In these early 
studies SF was traced by the associated Ho- emission (in HII re- 
gions) using interference filters, a technique widely used with 
this aim in the local universe since then. A considerable corpus 
of knowledge has been collected about the spatial distribution 
of SF in the local universe, characterizing its radial extent (e.g. 
Hodge [T969al Hodge & Kennicutt |1983l l, and studying how it re- 
lates to the morphological types of galaxies (e.g. Ho et al. |I997bl 
Shane |2002] |. the infall of new gas (e.g. Phookun et al. |1993l l, spi- 
ral pat terns ( e.g. Roberts I 19691 Cepa & Beckman 1990 Seigai- & 
James |2002] l, the influence of stellar bars (e.g. Phillips 1993. Ho 
et al. |I997bl l, interactions with other galaxies (e.g. Kennicutt et 
al. 1 1 9871 Moore et al. 1996), or with intra-group or intra-cluster 
media (e.g. Koopman & Kenney 2004|. 

In recent years there have also been some examples of ex- 
tensive studies on the radial distribution of SF in disc galaxies. 



with samples of dozens and even hundreds of galaxies. In The 
Ha Galaxy Survey, James et al. (2004) surveyed >300 galaxies 
as imaged in Ho-, and Shane (2002), basing his results on these 
data, studied the radial distribution (through several measures of 
the "concentration") of SF in interacting galaxies, in clusters and 
the field, and also as a function of the Hubble type. Other stud- 
ies have benefited from GALEX, such as Boissier et al. 120071 
They presented a study on the radial profiles in UV and IR for 
a sample of 43 nearby galaxies, and compared them to their dis- 
tributions of gas (CO and HI), studying the radial variation of 
SFR and dust-absorption in the UV. In another interesting ex- 
ample, Munoz-Mateos et al. 2007 studied the radial profiles of 
specific SFR, sSFR, for a sample of 161 nearby spiral galaxies. 
They found a large dispersion in the slope of these profiles with 
a slightly positive mean value, which they interpreted as proof 
of a moderate inside-out disc formation. Also based on deep 
GALEX observations has been the recent discovery of extended 
UV discs (XUV discs, Thilker et al. l2005l Gil de Paz 2005) . This 
UV emission, which is due to SF, and found far beyond the ex- 
tension of the optical disc, stands as a piece of evidence which 
is not easy to fit into the puzzle of disc formation. 

In this work we present a study of the radial distribution 
of rest-frame Near-UV flux in a sample of 260 disc galaxies, 
extending over the redshift range 0<z<l. The Near-UV flux 
(1500</1<2800A) is dominated by the contribution from young 
stars, with ages < lOOMyr (Kennicutt , 199 8), and so the flux in 
these wavelengths could be taken, to a fair approximation, as a 
tracer of on-going SF. Nonetheless, dust attenuation is signifi- 
cant at these wavelengths, and so, the observed distribution of 
NUV^ux does in fact reflect the distribution of SF, as modi- 
fied by extinction. In this work we explore the dust problem, 
although we are limited by the complexity of the issue, and lack 
of relevant data, particularly in the rest-frame infrared. We per- 
form some simple tests, based on known radial profiles of dust 
attenuation, to quantify the impact of dust extinction on the re- 
trieved results. Our aim is to probe the evolution in the spatial 
distribution of SF across stellar discs, and understand the results 
in the context of disc formation and evolution. 

This paper is organized as follows. In section |2] we describe 
our sample of galaxies, the data that we use and how they are 
handled for our study. In section [3] we present results on the ra- 
dial distribution of rest-frame Near-UV flux in these galaxies, 
and compare it to the distribution of flux in rest-frame B -band. 
Also in this section we analyse our results, accounting for the 
observational effects. In section|4]we discuss the results, relating 
them to the distribution of star formation, and put them in con- 
text, by relating them to other observational studies, and also to 
current models on disc formation. Finally, in section|5]we sum- 
marize our conclusions on the presented analysis. 

Throughout this work, we assume a flat A-dominated cos- 
mology (Q.M =0.3, Qa = 0.7, and Hq = 70 km s'^ Mpc^'^). 



2. Data & Selection of Samples 

The observational basis for this work is taken from several public 
databases which provide photometric and imaging data in rest- 
frame NUV and other bands for large samples of galaxies in the 
local universe and at intermediate redshifts (z < 1). Here we ex- 
plain the criteria for selecting the objects under study in different 
redshift bins, and provide information on the resulting samples 
and the data we use for the analysis. 



AzzoUini, Beckman & Trujillo: Radial Distribution of NUVFlux in disc Galaxies 



3 



2. 1 . Local Sample 

To build our Near UV (NUV) sample in the local universe we 
have chosen as parent sample The Galex Ultraviolet Atlas of 
Nearby Galaxies (ANG hereafter, Gil de Paz et al. 2007 ). This 
is composed of 1034 nearby galaxies, and includes objects from 
the GALEX Nearby Galaxies Survey (NGS ; Gil de Paz et al. 
120041 Bianchi et al. l2003all2003bl l. plus galaxies serendipitously 
found in the NGS fields or in fields with similar or greater depth, 
obtained as part of other GALEX imaging surveys, and which 
have optical diameters at the fig = 25 mag/arcsec^ isophote 
larger than 1 arc minute. 

To test for evolutionary changes, it is important to know to 
what extent this parent sample is representative of the galax- 
ies found in the nearby universe. With this aim, we refer to a 
comparison presented in Gil de Paz et al. 120071 of their sam- 
ple with the Nearby Field Galaxy Survey by Jansen et al. 120001 
(NFGS). This survey (listing 193 galaxies) was designed to be 
representative of the nearby population of galaxies, as portrayed 
in the magnitude-limited sample of the CfA Survey (Huchra et 
al. [T983] l. The GALEX -ANG and the NFGS are sampling the 
same volume of the universe and represent a similar population 
of galaxies, according to their distributions in redshift and ap- 
parent magnitude. The distributions of morphological types are 
also similar in these samples. There is, though, a small differ- 
ence in the luminosity distribution between samples. The most 
luminous spirals (-21<Mi, <-20) are in moderate excess in the 
ANG with respect to the NFGS , while the low-luminosity spi- 
rals {-l9<Mb <-18) are slightly under-represented. 

We define our sample in the local universe by applying the 
following criteria: a) it is restricted to disc galaxies; b) the ob- 
jects must be of moderate inclination, to facilitate the extraction 
of useful radial brightness profiles, and minimize the complica- 
tions in the analysis due to dust absorption; c) the objects must 
be bright enough for their counterparts at higher redshifts to be 
accessible for study with the available data; and d) we also want 
to have imaging data of the objects at longer wavelengths, within 
the optical range. With these criteria, we select the objects whose 
parameters are within these ranges: a) de Vaucouleurs morpho- 
logical type, T, within 0<r<10 (Sa to Sm) ; b) axial-ratio ^>0.5 
(inclination < 60deg) ; c) absolute magnitude in the B-band 
Mh <-19.5 mag (luminous galaxies above the completeness limit 
in the highest redshift bin, 0.8<z<l.l);d) not classified as part 
of a "multiple" system in HyperLeda; and e) within the footprint 
of the Data Release 6 of the Sloan Digital Sky Surve;^ (SDSS , 
Adelman-McCarthy et al. 2008) . The morphological types, axial 
ratios and Mb are taken from the HyperLedc0 database (Paturel 
et al.|20p). 

These criteria constrain our local sample to a population of 
98 nearby, almost face-on disc galaxies, extending to ~200 Mpc 
(z ~ 0.05). In addition, we further restrict the sample to galaxies 
which lie at a maximum distance of 60 Mpc. That is to avoid 
a bias in the selection that favours the largest galaxies, due to 
the criteria applied to build the ANG , our parent sample. In Fig. 
[T] we present the distribution of physical sizes (kpc) for galax- 
ies in the Local sample, as a function of distance (Mpc). These 
sizes are taken from HyperLeda, and correspond to the radius of 
the isophote at 25 mag /arcsec^ in the B-band. The dotted line 
marks the projected size of half an arcmin at the given distance. 
As mentioned above, the ANG contains objects serendipitously 
found in the NGS and other GALEX images which have a di- 



' http://www.sdss.org 
^ http://leda.univ-lyonl.fr 



ameter (at the 25 mag /arcsec^ isophote) equal to or larger than 
1 arcmin. Given that roughly 70% of the galaxies in the parent 
sample were included in this way, this means a selection effect 
by which there is a lower envelope in the distribution of sizes 
against distance, as observed in this figure (in Section[33]we jus- 
tify quantitatively the election of the imposed limit on distance, 
when the eff'ective radii of the galaxies are discussed). After this 
restriction in distance (D <60 Mpc), our local sample is com- 
posed of 33 galaxies. 

From this set we have discarded 2 more galaxies: NGC1042 
because of sky gradients in the SDSS -g' image; and NGC2782 
as suspect of undergoing a merger. This final sample of 3 1 ob- 
jects will be referred to hereafter as the "Local" sample. 

The dataset for exploring the radial distribution of flux in 
our Local sample is composed of images in 2 bands, which we 
list with their corresponding effective wavelengths (/If//) and 
full widths to half maximum of the transmission curves (A/l): 
GALEX -NUV {Aeff = 2250A, = lOOOA) and SDSS- 
g'(Aeff = 4770A, AA = 1370A). The GALEXimages have a 
PSF which is slightly broader for brighter sources, and vari- 
able throughout the field of view (FOV), which is circular and 
1.2deg in diameter. For objects within the central 0.5deg of the 
GALEX FOM, the full width at half maximum (FWHM) of the 
Point Spread Function (PSF) is in the range 5.0"-5.5" in the 
NUV hand. The projected pixel scale on the sky is 1.5" per side. 
The SDSS -g' images have a median value of the FWHM of 
1.5", and an angular scale of 0.396"/pixel. 

The ANG images have a characteristic exposure time of one 
orbit (~1700 sec), while for the SDSS images the exposure 
time is fixed at 54 sec. Together with other relevant parameters, 
in Table [1] are listed Icr fluctuations in surface-brightness of the 
background in an aperture with an area of 1 arcsec^, in all bands 
used in this work io-„rcsec-, 6th column; all magnitudes there and 
throughout this work are in the AB system). We also list the 
Icr fluctuations of the background in the typical area used to 
estimate the background, an elliptical annulus with semi-major 
axis equal to 1.5 times the isophotal radius of the object in the 
B/tf band, and width of 1 pixel {a-i,ack, 7th column). This isopho- 
tal level to which we refer is, of course, that used to detect the 
objects. As we can see in this table, there are substantial diff'er- 
ences between the noise levels amongst bands. How we account 
for these differences in our analysis is explained in detail in sec- 
tion |3] 

We could have benefited from other bands also available 
both at shorter and longer wavelengths to extend the analysis 
on the radial distribution of flux for galaxies in the Local sam- 
ple. To extend to the "blue" side, which would be interesting, 
for example, to estimate dust absorptions in the UV, we could 
use the GALEX -FUV images {A^ff - 1520A), that are avail- 
able for most of the galaxies in the ANG surveys. On the "red" 
side, which could provide information on the distribution of yet 
older stars than those traced with rest-frame B -band, r' , /' and 
z' SDSS bands (/i?// > 6000A), are also available. The reason 
not to use them is that we are limited by the available bands at 
z~l from GOODS, if we want to make a meaningful compari- 
son of the radial distributions of flux of the objects. This is, we 
make this restriction in order to compare distributions of fluxes 
in similar rest-frame wavelengths at all redshifts. Nonetheless, 
we use all bands listed in Table [1] including all iheSDSS bands 
for the Local samples, to make estimates of B -band luminosities 
and stellar masses of the galaxies, by means of SED modelling, 
as described in section |3] 
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Fig. 1. Radii of the galaxies at the 25 mag / arc sec~ isophote in 
the B -band (R25) against distance, for objects in the Local sam- 
ple. The dotted line corresponds to the length subtended by half 
an arcmin as a function of distance, which is a constraint put on 
the selection of many of the objects in the ANG survey. The ver- 
tical dashed line marks a distance of 60 Mpc, which is a further 
restriction we set on the sample to minimize the bias in apparent 
size inherited from the ANG . 



2.2. Intermediate Redshift Sample 

We survey galaxies which extend up to z~l, which corresponds 
to a look back time (Ibt) of ~8 Gyr, according to the most 
favoured A-CDM cosmological model. In order to have an imag- 
ing database which provides rest-frame A^J/V coverage, and with 
usable spatial resolution and sensitivity levels, we chose the 
HS T -ACS imaging of the GOODS -Soutlfl field (Giavalisco et 
al.|2004). This field, which has been intensively observed using 
several space and ground observatories (e.g. HST , Chandra, 
XMM - Newton, S pitzer, VLT, etc.), subtends an area on the 
sky of roughly 170 square arcmin. The HST -ACS imaging set 
consist of images in the F435W, F606W, F115W and FE50LP 
HS T pass-bands, hereafter referred to as B435 , Veoe , '775 and 
Z850 , as is common in the literature, and for which total expo- 
sure times per pixel in the final images were of approximately 
7200, 3040, 3040 and 6280 seconds, respectively. The images 
have an angular scale of 0.03"/pixel and the FWHM of the PSF 
in zg5o measures 0.09". 

We have benefited from the fact that the GOODS -South 
field is one of the most thoroughly observed patches in the sky, 
and thus we rely amply on the work of other authors to define our 
sample at intermediate redshift and gather valuable data on our 
targets. First, the GOODS -South field is contained in the Galaxy 
Evolution from Morphologies and SEDs imaging Survey (GEMS 
; Rix et al. l2004l) . In the redshift range 0. 1 <z< 1 , GEMS provides 
morphologies and structural parameters for nearly 10,000 galax- 
ies, obtained using //ST images (Barden et al. 2005; Mcintosh 
et al. l2005]) . Second, for many of these objects there exist pho- 
tometric redshift estimates and Spectral Energy Distributions 
(SEDs) from COMBO-17 {Classifying Objects by Medium- 
Band Observations in 17 filters; Wolf et al. 2001, 2003). The 
COMBO-17 team made this data publicly available through a 
catalogue with precise redshifts (with errors &/(l-i-z)~0.02) for 
approximately 9000 galaxies down to m«<24 mag (Wolf et al. 
[2004), which we use. In the same data release were included 
rest-frame absolute magnitudes and colours (accurate to ~ 0.1 
mag). Finally, we have also used the stellar mass estimates pub- 



lished in Barden et al. 120051 which are taken from Borch (I2004I I. 
and are deduced from the COMBO-17 photometric data. 

Barden et al. (2005) conducted a morphological analysis of 
the galaxies in the GEMS field by fitting Sersic r^l" (Sersic fT968] l 
profiles to the surface brightness distributions. Ravindranath et 
al. (|2004) showed that using the Sersic index n as the criterion, it 
is feasible to distinguish between late and early-type galaxies at 
intermediate redshifts. Late-types (Sab-Sdm) are defined as hav- 
ing n <2-2.5. Moreover, the morphological analysis conducted 
by Barden et al. provides information about the inclination of 
the galaxies, through the axial ratio of the isophotes, q. This is a 
parameter to take into account, since we want to study the radial 
distribution of the rest-frame NUV flux. This is particularly sen- 
sitive to dust absorption, and thus to the inclination of the discs, 
as already mentioned in the description of the selection of the 
Local sample. 

Our intermediate redshift sample, to which we refer through- 
out this paper as the High-z sample, is constructed from the col- 
lection of data just described as follows. First, we selected ob- 
jects from the Barden et al. sample within the following ranges 
of parameters: Sersic index n <2.5 to isolate disc-dominated 
galaxies (Barden et al. 120051 Shen et al. 120031 Ravindranath et 
al. |2004); axial ratio q >0.5 to select objects with inclination 
<60deg ; andM^ <-19.5 mag, as for the Local sample. Moreover, 
only objects within 0.5<z<l.l were selected in order to have 
spectral coverage in the rest-frame At/Vband with the avail- 
able images from HST -ACS. Then, the resulting sample was 
matched to a photometric catalogue extracted by ourselves from 
the GOODS -South HS T -ACS data {GOODS data hereafter), 
using SExtracto:0 (Bertin & Arnouts 1996). The resulting sam- 
ple contains 265 objects. Also in this case it was necessary to dis- 
card a number of objects because of diff'erent problems that make 
them ill conditioned for our analysis: a) proximity to a bright 
source, such as an star or galaxy, b) image artefacts, a condition 
that is more relevant at the borders of the images, c) the object is 
apparently undergoing a major merger event. After these rejec- 
tions, the High-z sample is composed of 239 disc galaxies. 

Furthermore, we divide the High-z sample into two subsam- 
ples, according to redshift: a) "mid-z" subsample, with 134 ob- 
jects within 0.5< z <0.8, and b) "far-z" subsample, with 105 
galaxies within 0.8< z < 1.1. In the first redshift range, the 
GOODS band that best traces the rest-frame B -band is ijj^ , 
while in the second range, it is zsso ■ We introduce this segre- 
gation to allow for a more accurate comparison of the radial dis- 
tribution of At/yflux relative to rest-frame B-band flux, as the 
changes in observing band due to redshift are less significant in 
this way. 

In Fig.|2]we show the distributions of Z?-band luminosities 
{Mh ), and stellar masses (M* ) for objects in the Local and High- 
z samples. The luminosities are quite similar amongst Local and 
High-z galaxies, with median values of -20.4 mag and -20.6 mag 
respectively. The distributions of masses are not so similar, with 
the Local sample extending somewhat to more massive galax- 
ies. The Local sample has a median value of stellar mass, , 
of log(M* )=10.1, and the High-z sample of 9.6. These diff'er- 
ences have relevant implications for the analysis strategy, as we 
are interested in comparing the radial distribution of NUV flux 
amongst galaxies which are indeed comparable, so as not to con- 
found selection effects with genuine evolution. We explain this 
strategy extensively in section [3] The luminosities and stellar 
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This catalogue was obtained by detecting sources in the zgso band 
which had at least 16 contiguous pixels (0.014 arcsec^) at an isophotal 
level of 0.6 cr j(.,,/pixel (25.4 mag I arc sec- ) or higher. 
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Table 1. Relevant Parameters of the Imaging Dataset. Columns: 1) band; 2) effective wavelength of the transmission curve; 3) 
full width at half maximum of the transmission curve; 4) typical exposure time; 5) projected size of the pixel in the sky; 6) 1 cr 
fluctuations of the background in an area of 1 arcsec^; 7) 1 cr fluctuations of the background in the typical area used to estimate the 
background (an elliptical annulus with semi-major axis equal to \.5Rp (Petrosian radius in the Bs/rband), and width of 1 pixel); 8) 
FWHM of t\\s PS F . 

Notes: a) NUVrf at z ~0; b) Brf at z ~0 ; c) NUVrf within 0.5< z <0.8 (mid-z) and 0.8< z <1.1 (far-z); d) Brf within 0.5< z <0.8 ; 
e) Brf within 0.8<z<l.l. 
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Fig. 2. In the upper panels are shown the distributions of 
Luminosities (M^ , left), and stellar masses (M* , right) for 
galaxies in the Local sample. In the lower row are the corre- 
sponding distributions for galaxies in the High-z sample: dotted 
line for the mid-z subsample, dashed for the far-z subsample, 
and continuous for the whole High-z sample. 
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Fig. 3. Left: Morphological types (de Vaucouleurs) of galaxies 
in the Local sample. Right: Sersic index of the galaxies in the 
High-z sample (continuous line). Dotted line: mid-z subsample 
(0.5< z <0.8); dashed line: far-z subsample (0.8< z <L1 ). 



masses shown in this figure for both samples, and used through- 
out this work, are based on our own photometry of the objects, 
and obtained through model fitting of the SEDs with HyperZ, as 
described in section [3] 

Another characteristic of the galaxies which is relevant to 
our study is the morphological type of the galaxies, as different 
trends have been reported in the spatial distribution of SF for 
different Hubble types of galaxies (e.g. Shane 2002). More im- 
portantly, one of the aims of this study is to provide new data 
to help understand how the location of newborn stars is related 
to the galaxy morphology. So, it is imperative to select galaxies 
with similar morphologies in the different redshift bins. In Fig.|3] 
we present the distributions of morphological types for objects in 
the Local sample (left; taken from HyperLeda) and of Sersic in- 
dex for objects in the High-z sample (right; Barden et al. 120051 1. 
The morphological type is related to the Sersic indexes (Barden 
et al. 2005, Shen et al. 2003, Ravindranath et al.'200T), and un- 
der this premise we can compare the morphologies of the two 
samples with the available data. Ravindranath et al. j2004l) ex- 
plored the relation between the Sersic index and morphological 
type of galaxies when the effect of redshift on their observation 
is simulated. They retrieved the Sersic indexes for local galax- 
ies artificially redshifted to z~ 1 (as they would be observed in 
the GOODS-S field). In Fig. 1 in their paper it can be seen how 
galaxies with 0.5 < n <2, the range in which most of the objects 
in our High-z sample are found, have preferably morphologi- 
cal types, T, around 5, fairly similar to the most frequent mor- 
phological type in the Local sample. Thus, the morphological 
ranges are quite similar amongst low and intermediate redshift 
samples, and the most frequent morphological types are in the 
range 2<T<6, i.e. Sab-Scd. 



2.3. Preparation of the Images for Analysis 

Due to the heterogeneity of our data, we need to perform several 
operations on the images before the analysis is done. Here we 
describe these operations, first for the Local sample and then for 
the High-z sample. 

The GALEX -NUV images from the ANG are already sky 
subtracted, calibrated and registered to a World Coordinate 
System (WCS), and the objects of interest are fully included, 
as it is an atlas centred on targets, not on coordinates. In con- 
trast, the same galaxies in the SDSS are in many cases not in 
the field centre, and may even be split amongst two or more 
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"fields" (using the SDSS nomenclature for the released indi- 
vidual images), depending on size and position, and the im- 
ages are not sky subtracted (though SDSS provides estimates of 
the background values). We thus have to produce our own sky- 
subtracted mosaics of the galaxies from the S DSS images. The 
background is estimated from the median of counts in regions 
where there are no detected objects, in an iterative procedure, 
for more accuracy. The mosaics produced are WCS-registered 
to the GAL/iX images, but with the angular scale of the S DSS , 
0.396"/pixel. 

One problem that arises is that in these images there may be 
contributions by stars from our own galaxy, and also field galax- 
ies, whose presence is spurious to our aims. We remove these 
intruding objects in two steps. In a first step, we create a mask 
of stars from the USNO-B star catalogue (Monet et al. 2003). 
On the position of each star, and in a circular area with diameter 
equal to a certain factor,^*, times the full width at half maxi- 
mum (FWHM) of the PS Fin the g'band, we have flagged the 
pixels as "Not a Number", NaN, thus being discarded from the 
computations. The aforementioned factor was fixed at =2 as 
a best solution, after some trial and error. This mask is visually 
inspected and manually edited to ensure that in a sensitive region 
on and around the object of interest, all stars, and only the stars, 
are eff'ectively being masked. This is necessary as the USNO-B 
catalogue is not always complete in accounting for all detected 
stars at the depth of the images, and besides, some bright blue re- 
gions in the galaxies (zones with recent or ongoing intense star 
formation) are in cases wrongly classified as stars in the cata- 
logue, and we do not want to blank out these zones, given the 
nature of our study. In a second step, we manually create a sec- 
ond mask, to blank out those regions of the images (i.e. set them 
as NaN, thus being ignored altogether) which are occupied by 
extended field galaxies or bright stars (whose bright wings of 
the PSF are not effectively masked in the first step). This mask- 
ing does not prevent some galaxies from having to be discarded 
as not suitable for the study, as we mentioned above, when the 
impact of intrusive objects is too severe, or there are other prob- 
lems with the quality of the images. 

With the High-z sample the actions previous to the analysis 
are somewhat less involved. In first place, the GOODS images 
are already registered to a common WCS, photometrically cali- 
brated and sky subtracted. So we start dividing the images into 
"stamps" with a typical size of 5 arcseconds on a side, to speed 
up the computational tasks. Also here too there are spurious ob- 
jects in the neighbourhood of the galaxies under study. But in 
this case, given the high galactic latitude location of the field, and 
given that it subtends a relatively small area on the sky, there are 
few stars of our galaxy in the image stamps of the objects. The 
presence of field galaxies is more important, in number, because 
of the depth of the images. In this case these intrusive objects 
are masked by blanking out those areas of the images which are 
classified by SExtractor as belonging to any object other than the 
detected target. We perform a visual inspection of each stamp af- 
ter this masking to ensure that it looks reasonable, and discard 
those objects for which it does not work properly, as is the case 
when there is contact (apparent or physical) between the intrud- 
ing objects and the target. This "SExtractor masking" procedure 
is also performed for the local objects, but in many cases it turns 
out to be inefficient in eliminating the stars, especially when they 
are superposed on the object of interest. That is why for objects 
in the nearby sample we need to apply, in addition, the more 
involved procedure described above. 



2.4. Differences in relevant observational parameters 
amongst the imaging datasets 

Another problem that arises when one aims at comparing the 
radial distribution of luminous flux at a given wavelength of 
galaxies at different redshifts, and observed through different in- 
struments, is that, in general, the images used for this task trace 
those distributions in different ways, which may bias the study. 
The most important parameters which may differ amongst the 
images are: a) angular scale (y, given in arcseconds/pixel); b) 
physical scale, defined as the proper length that the angular scale 
subtends at the distance of the object (F, given in kpc/pixel) c) 
PSF properties (most importantly its FWHM ) ; d) photometric 
depth ; and e) relative band-shifting due to differences in red- 
shift. Variations in these parameters amongst images makes the 
direct comparison of the flux distributions as registered in them 
misleading, or at least, more difficult to interpret. In this work 
we take into account these issues when interpreting our results, 
as we explain in this subsection. 

2.4.1. Resolution 

The objects in the Local sample extend in distance between 
4 and 60 Mpc. With the GALEX angular scale (y=1.5"), this 
means a physical scale F varying in the range 0.03 < F <0.44 
kpc/pixel, while with the 5 5 pixel (0.396") the same pa- 
rameter ranges between 0.008 and 0.12 kpc/pixel. On the other 
hand, the small angular scale of the HST -ACS images (y = 
0.03"/pix after "drizzling", i.e. with high spatial resolution), 
translates, within the redshift range 0.5<z<l.l, into a physical 
scale 0. 18< F <0.25 kpc/pixel. So, the HS T images in fact pro- 
vide a better resolution per pixel than GALEX for the farthest 
objects in the respective samples, though for the nearest galax- 
ies, the SDSS and GAL/sX images provide significantly higher 
resolutions (roughly by an order of magnitude). 

With regard to the PSF of the images, we assume it to be 
Gaussian, in all cases to a good approximation, with the only 
parameter necessary to describe it being the FWHM . This pa- 
rameter also varies over a wide range amongst the datasets. For 
GALEX -NUV images the FWHM is ~5", and this is the value 
we adopt for these images. For the S DSS images this parame- 
ter is, in general, different for every object, because the images 
were taken under varying seeing conditions. A guide value in the 
g' band would be ~ 1 .5". For this dataset we use, in each case, the 
value of the FWHM in the Field image where the centre of the 
object lies, in every band, as reported in the corresponding ts- 
Field file (again using the SDSS nomenclature). For the High-z 
sample, we assume a FWHM of 0.09" for all bands, as the minor 
differences amongst them are not significant to our study. 

As with the size of the pixel, there are also significant vari- 
ations in the projected dimension of the FWHM , measured in 
kpc, at the distance of the galaxies. For example, for an object 
at 60 Mpc the GALEX FWHM of 5" subtends -1.45 kpc, while 
the 0.09" FWHM of the GOODS images ti-anslates into a better 
resolution of F ~0.7 kpc at z=l. 

In this paper we present an analysis of effective radii (i.e. the 
major axis of the -assumed elliptical- isophote which encloses 
a 50% of the total flux, R^q ), and of "Concentration", in whose 
expression appear Roo and /Jgo , which enclose 20% and 80% of 
the NUV flux respectively. In order to show that we can measure 
the above quantities accurately, we compare these radii with the 
FWHMs. In Fig. |4]we present the ratio of these radii (referred 
to A^f/y band) to the FWHM of the PSF in NUV , for different 
redshift bins. We use this code to differentiate the samples: blue- 
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continuous line for the Local sample (z~0), green-dotted for the 
mid-z subsample (0.5< z <0.8 ) and red-dashed for the far-z sub- 
sample (0.8< z < 1 . 1 ), and this same code is used throughout this 
work in histograms of properties of these samples. In the left 
panel we see the aforementioned ratio for R20 , which would be 
the radius most affected by differences in spatial scale as the red- 
shift changes. The worst cases are for galaxies in the mid-z and 
far-z subsamples, the minimum values of the ratio 2-/? 20 /FWHM 
being 0.9 and 1.7 (for two objects at z =0.57 andz =1.07 respec- 
tively). The median values are between 5 and 9 (i.e. large enough 
to not be affected by the PSF), with the objects in the High-z 
sample having comparatively the worst resolutions. In the other 
two panels the median values of this ratio range between 10 and 
15 for /?5() and between 17 and 23 for R^o ■ So, the dimensions 
which are of interest to us, these radii, are sufficiently larger than 
the FWHMs, though there is a non negligible difference between 
the effective resolution for the Local and High-z samples. 

The figures just reported could raise doubts about whether 
our analysis is somehow biased due to differences in spatial res- 
olution with redshift and datasets. To double check this issue 
we also study what happens when the images of the objects are 
degraded to the same spatial resolution. This common resolu- 
tion is given by the lengths, in kpc, that the pixel and FWHM 
of the GALEX images subtend at 60 Mpc, the maximum dis- 
tance for galaxies in the Local sample. At that distance, the 
angular scale of 1.5"/pixel gives a physical scale of F =0.44 
kpc, and the FWHM of 5" subtends 1.45 kpc (we dub this 
magnitude FWHMkpc )• By means of interpolation and convo- 
lution, we produce images in rest-frame B -band, Brf , and rest- 
frame NUV -band, NUVrf , of all galaxies under study at that 
resolution (after rounding off to F =0.5 kpc, FWHMkpc - ^-^ 
kpc). How this is done is explained in the Appendix. In Fig. 
|5] we can see examples of this transformation for two galax- 
ies, at low redshift (Messier 95, D=12 Mpc ; above) and in- 
termediate redshift (GEMSz033236.03m274423.8, z=0.95 ; be- 
low), presenting their "original" resolution (left), and degraded 
to the common resolution (F =0.5 kpc, FWHMkpc =1.5; right). 
Throughout this work we comment on how the use of images 
with a common resolution affects the reported results. 

2.4.2. Band-shifting 

Another key issue to take into account in our study is that as 
we observe objects in a range of redshifts, the observing bands 
trace different parts of the rest-frame SEDs of the objects. This 
must be accounted for if meaningful results on the distribution 
of the underlying stellar populations are to be obtained. We have 
decided on a strategy such that the selection of observing bands 
in which to extract the different measures depends on the red- 
shift "bin" in which the object lies. This is a relevant issue for 
the High-z sample, but the observing bands in the Local sam- 
ple must be taken into account as well to allow for comparison. 
As we describe further on, several parameters of the galaxies are 
measured in a band selected to trace the rest-frame B -band, at all 
redshifts. Thus we refer to it as the band. In a broad sense, 
we use this band to characterize the distribution of "older" stars 
in the galaxies, to complement and contextualize the descrip- 
tion of those newborn/very young stars, traced by the rest-frame 
NUV images. For the Local sample, this Brf band is the g' band. 
In the High-z sample, it depends on the redshift bin. For the mid- 
z subsample (0.5< z <0.8) it is the /775 band,and for the far-z 
subsample (0.8< z <1.1 ) it is the Ziso ■ In a similar way we also 
define a NUVrf band to refer to the band that we use to trace 
the rest-frame A^f/V flux in each redshift bin. For objects in the 



Local sample, this corresponds to the G ALEX -NUV h&n^^ itself, 
and for those in the High-z sample it is the ACS-B435 band. 

2.4.3. Pliotometric Deptli 

We also must discuss the differences in depth between the 
datasets, which are significant. To better understand these 
differences we will give some illustrative numbers. The 
Icr/pix noise levels, and angular scales, y, for the GALEX , 
SDSS and GOODS data, in the A^f/Vs/r band, are respec- 
tively 28 magjarcsec^ - 1.5", 24 magjarcsec^ - 0.396" and 25 
magjarcsec^ - 0.03". Then we assume characteristic distances 
for each sample: 30 Mpc for galaxies in the Local sample and 
z =0.8 for those in the High-z sample. If we compute how many 
pixels are needed on the given images, at the corresponding dis- 
tances, to cover an area of 1 kpc^ we get these figures: GALEX , 
21 pixels; SDSS , 301; and GOODS , 20. If we take the given 
noise levels, and take into account that when we average the 
intensities over N pixels the signal-to-noise ratio increases by 
the V^, then the corresponding Icr/kpc^ values would be 29.7 
mag/kpc^ (GALEX, Local), 27.1 mag/kpc^ (SDSS, Local) and 
26.6 mag/kpc^ (GOODS, High-z). But at z =0.8 we have 2.6 
magnitudes of cosmological dimming, and so, the actual value 
to consider for GOODS-High-z would be 24 mag/kpc^. So, on 
equivalent areas, GALEX goes ~10 times "deeper" than SDSS, 
and SDSS goes <20 times deeper than GOODS (cosmological 
dimming included). In section |4] we discuss how this may affect 
our results. 

3. Data Analysis & Results 

3.1. Extraction of Radial Surface Briglitness Profiles 

The main goal of this work is to trace the evolution of the radial 
distribution of SF in disc galaxies, using as proxy for this dis- 
tribution that of the rest-frame NUV^nx. To accomplish this, 
we use the analysis of surface brightness {p - r) profiles in 
different bands. These p - r profiles were extracted using pho- 
tometry on quasi-isophotal elliptical apertures. The intensities 
were estimated as the median of the flux in the area between 
elliptical apertures of increasing semi-major axis length (here- 
after we term these lengths, loosely, as radii). The ellipticity and 
position angle of the apertures are fixed to those retrieved by 
SExtractor for the whole distribution of light of the object in the 
Brf band. The pixels that go into this calculation are those in the 
Brf isophotal area of the object (see below for an explanation 
of this). The centre of the apertures is also fixed. For the Local 
sample, the centre is that reported in HyperLeda for the object, 
which after inspection seems adequate in all cases. For the ob- 
jects in the High-z sample the process is more involved, because 
these objects show, in many cases, appreciable asymmetries. The 
procedure to pin-point the centre of the apertures is this: in a first 
iteration, the first moments (in "x" and "y") of the brightness 
distribution of flux of the object in the zgsoband are employed. 
Z850 is the reddest band available for the High-z sample, and thus 
traces best the distribution of old stellar populations, and thus, 
stellar mass. After visual inspection of the resulting profile and 
the image of the object, the centre was refined, if needed, with 
help of the task "imexam" from IRAF0, to match what was vi- 
sually estimated, in each case, as the dynamical centre of the 
object. In a regular disc galaxy, as are those under study, this 
centre coincides with the central bulge or nucleus. 

' http://iraf.noao.edu 
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Fig. 4. Quantifying the resolution of our images: ratio of the diameter of the isophote which encloses a certain fraction of the total 
flux of the galaxy in NUV (20%, 50% and 80% from left to right) to the FWHM of the PS F , depending on the redshift range. The 
samples and subsamples are coded as: blue-continuous for the Local Sample, green-dotted for the mid-z subsample (0.5< z <0.8 ) 
and red-dashed for the far-z subsample (0.8< z <1.1 ). For objects of the mid-z and far-z subsamples, the FWHM is that in the 
£435 band. The vertical lines mark the median value of each distribution, and this applies to all other histograms shown in this 
article. 






Fig. 5. Examples of projections of images to a common resolution. The galaxy portrayed on the top line is Messier 95 {D -12 
Mpc) in GALE X -NUV hand, at its intrinsic resolution on the left (y =1.5"/pixel), and with degraded resolution (pixel =0.5 kpc, 
FWHM =1.5 kpc), on the right. On the bottom line we show the analogous sequence for GEMSz033236.03m274423.8 (z=0.95), 
in band GOODS -B435 . The red bar marks a length of 5kpc at the distance of the galaxies. 



The radius of the annular apertures was linearly increased 
in steps of 1 pixel (the corresponding angular scale varies de- 
pending on the dataset; see Table[T]for reference), up to a radius 
which is 50% larger than the radius of a circle with the isophotal 
area of the object in the Z?R/rband. The isophotal area is that 
covered by the set of connected pixels with intensities above 
the detection threshold (Icr/pixel) which constitute a detection 
(in our case: 24.7 mag/arcsec^ in g' , 25.2 mag I arc secern i-j-j^ , 
and 24.8 mag/arcsec^ in zs5Q , for objects in the respective sub- 
samples). For objects with a more or less regular morphology, 
as are those selected, the median intensities in these annulii are 
a good approximation to isophotal intensities. The error in the 
intensity, 61, is given by the cr of the distribution of fluxes in- 
side the annulus, divided by the square root of the correspond- 
ing pixel area. The intensities ("I") are transformed to surface 
brightnesses (yu), expressed in mag / arc sec^ in the AB system, 
using the corresponding magnitude zero point and angular scale, 
y, through the expression fi - zero - 2.5-log{l/y^). For the errors 
in ju, 6fi, the formula used is d/u - 2.5-log( 1 + 61/T). We produce 
surface brightness profiles of the objects in all available bands, 
according to the scheme described. For this task we use several 



scripts written by us in Python languag^ which "glue" together 
different tasks of IRAF, SExtractor and DSSQ 

The previously described intensity profiles are used in the 
analysis of stacked profiles (sect. 13. 4t , and in the estimate of 
centre-to-border differences in surface brightness (sect. 13. 7| ). To 
measure effective radii, and other magnitudes which rely on the 
amount of flux contained in different radii (e.g. the concentra- 
tion parameter, C42 ), we recur to the growth curves of the ob- 
jects. These are computed by integrating the flux in the elliptical 
apertures previously described. The corresponding radii are the 
semi-major axes of the elliptical apertures. 

In Fig.|6]we show several examples of objects in the Local 
and High-z samples, as seen in rest-frame visible bands (as 
"rgb composites", on the left), in rest-frame NUV (centre), 
and also present their surface brightness profiles in all avail- 
able bands (left). In the profiles we see clear truncations, also 
termed Type II profiles (Freeman 119701 Erwin et al. I2005I I. 
in NGC3319 (D =15 Mpc, first from above) and NGC3359 
(D =18 Mpc, second). One object shows a Type II profile 

* http://www.python.org 

' http://heq-www.harvard.edu/RD/ds9/ 



AzzoUini, Beckman & Trujillo: Radial Distribution of NUVFlux in disc Galaxies 



9 



(GEMSz033237.54m274838.7, z =0.67, fourth), though this 
"downbending" seems not to be a genuine disc truncation, 
as it happens well inside the disc. Another case also shows 
hints of truncation (GEMSz033251.15m274756.0, z =0.54, 
third), though not so clearly, as it happens in a region of 
the profile with low signal-to-noise ratio. The fifth object, 
GEMSz033248.47m275416.0 (z =0.78), shows a single expo- 
nential profile (Type I). As we can see, the NUV images show in 
all cases a patchy appearance, with varying degrees of similarity 
to the visible images. In the local objects it is easier to recognize 
structures such as spiral arms in the A^f/VR/r images, and also 
in GEMSz033237.54m274838.7, though in this case the spiral 
arms show some interesting asymmetry, perhaps the signature of 
a past interaction. The other two galaxies, taken from the High-z 
sample, show less defined morphologies, though the presence of 
a stellar disc is clear in all cases. 

As we have already stated, we use the A^t/V surface- 
brightness profiles and growth curves as proxies for the distribu- 
tion of the SF in galaxies. Nonetheless, the dust extinction plays 
a significant role in the shaping of these profiles that cannot be 
ignored, as the extinction curves usually adopt a power-law pro- 
file in which there is higher absorption at shorter wavelengths 
(A{A) oc A^^ to a first order approximation). A detailed treatment 
of this issue would require us to produce radial extinction curves 
for every object, which is by no means a straightforward task. 
To be done reliably it requires more data than we have (it would 
be desirable at least to extend the spectral coverage into the rest- 
frame near IR). Present state-of-the-art instrumentation (i.e. the 
absence of high-resolution IR facilities) prevents us from tack- 
ling the problem in such depth. So, in the next subsections we 
present results which do not include extinction corrections, and 
so our picture of the radial distribution of the SF and its evolution 
since z ~l should be taken with caution. 

Nonetheless, as a first and crude approach to solving this 
problem, in Sect. 13.91 we explore how the results might be af- 
fected by a radially resolved correction on dust extinction, fol- 
lowing the work by Boissier et al. (I2004I I. At the very least that 
exercise serves to visualize that the effect of dust cannot be ig- 
nored in this kind of studies, and also as an indication of the 
sense in which more reliable corrections would alter the pre- 
sented results. 

3.2. B-band Luminosities and Steiiar Masses of Gaiaxies. 

As we stated in section|2] we have estimates of the B -band lumi- 
nosities, Mh , and stellar masses, , of the galaxies in the High- 
z sample from the COMBO- 17 dataset. We also need similar es- 
timates for the Local sample, so as to compare, amongst redshift 
bins, the relation between different measures of the radial dis- 
tribution of NUV -finx and these physical properties. With this 
aim we have, first, performed aperture photometry of the galax- 
ies in the Local sample, to build their SEDs in NUV and the 5 
5Z)55 bands. Then, using the HyperZ code (Bolzonella et al. 
I2OOO1) . we fit these photometric data to model SEDs to obtain es- 
timates of the desired properties of the galaxies. We use the stel- 
lar population synthesis models implemented in HyperZ, which 
are based on synthetic spectra from the GISSEL98 ("Galaxy 
Isochrone Synthesis Spectral Evolution Library") spectral evo- 
lution library of Bruzual & Chariot il993i The Initial Mass 
Function (IMF) used is that of Miller & Scalo |1979l with lower 
and upper cut-offs in stellar mass at Mi =0.1 Mq and M2 - 125 
Mq . With respect to the Star Formation History (SFH) we have 
opted for Single Stellar Populations (SSP) of different ages 
(from 1 Myr to the age of the universe at the redshift of the ob- 



ject). We assume a solar metallicity (Z = ~ 0.02) and the 
SEDs are attenuated, using a free parameter. A,,, according to 
the law of reddening by inner dust absorption given in Calzetti 
et al. 120001 for starburst galaxies. We restrict this parameter to 
be within 0<Av<2 mag. The reddening produced by dust in our 
own galaxy is also introduced in terms of the excess in the B -V 
colour for each galaxy {E(B - V); estimates on this parameter are 
taken from HyperLeda), given as input. The redshift, which is 
normally the parameter needed when using HyperZ, is fixed, and 
given by the corresponding distance according to HyperLeda. 

The reliability of part of the analysis presented in this work 
depends on how consistent are our previous estimates of stellar 
masses and luminosities of galaxies with those from COMBO- 
17. In order to check this, we proceed as follows. We have mea- 
sured the Ml, and Mi, of the High-z galaxies with HyperZ, in 
an analogous way to that described above for the Local sam- 
ple. The only differences in the input data are the photometric 
values (here we use the 4 available GOODS bands), and the red- 
dening from the Milky Way, which is fixed to be E(B-V)=0.007 
mag (Schlegel et al. I1998I I. Then we compare the estimates of 
Mb and thus produced with those from COMBO- 17, and the 
results are shown in Fig. [T] The Mh estimates (left) show good 
agreement, our values being 0.2 dex brighter than those from 
COMBO-17, with relatively little dispersion (cr=0.3 mag). The 

estimates have a somewhat larger dispersion (cr=0.2 dex), 
and there is some bias, our mass estimates being lower than those 
from COMBO-17 by 0.14 dex (-30% less massive) (at M^ = 
10'" Mq). The main contributors to this bias are, presumably, 
three factors: a) the photometry; b) the choice of IMF; and c) 
the SFH of the model SEDs. The photometry from COMBO-17 
was obtained from ground observations (MPG/ESO 2.2m tele- 
scope on La Silla), and applying aperture corrections depend- 
ing on morphological parameters to retrieve the fluxes, while we 
have just integrated the fluxes in the, elliptical, quasi-isophotal 
aperture with semi-major axis equal to the rest-frame B-band 
Petrosian radius of the object. Second, we use the IMF from 
Miller & Sc alo ( fT979] l while COMBO-17 uses that from Kroupa 
et al. (I1993I I. We have tried other SFHs, apart from single bursts, 
such as a constant SFR or an exponentially decaying SFR, but 
the best results, in terms of scatter and bias when comparing with 
COMBO-17, were those obtained with the SSP models, which 
are thus our choice. 

To conclude this point, and with aim of minimizing the pos- 
sibility of having any bias when comparing results which depend 
on the stellar masses or luminosities of galaxies at different red- 
shifts, throughout this work we use our own "HyperZ" estimates 
of these quantities at all redshifts, rather than the values from 
COMBO-17 available for the High-z sample. This includes the 
results shown in Fig|2]on the distributions of stellar masses and 
luminosities for the different redshift samples, and already com- 
mented in section |2] 

3.3. Effective Radii (R,ff ) 

The most basic characterization of the radial distribution of the 
Near UV flux is the effective radius. This is defined as the radius 
which encloses half the total flux of an object, given a certain 
geometry of apertures (i.e. their centre, ellipticity and polar an- 
gle, assuming they are elliptical). We use the same apertures em- 
ployed to retrieve the radial surface brightness profiles described 
in section|2] In Fig[8]we show the distributions of effective radii 
in NUVrp for galaxies in different reshift bins. The same com- 
binations of line style and colour that were used in Fig. |4] ap- 
ply here to code the different redshift subsamples, and this code 
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Fig. 6. Examples of images and surface-brightness profiles of objects under study. From above: NGC33 19 (D -15 Mpc), NGC3359 
(D =18 Mpc), GEMSz033251.15m274756.0(z =0.54), GEMSz033237.54m274838.7 (z =0.67) and GEMSz033248.47m275416.0 
(z =0.78). From left to right we show an rgb-composite, the image in NUVrf and the observed surface brightness profiles in all 
available bands (cosmological dimming has not been accounted for in these figures). The bands used as the 'r' 'g' and 'b' channels 
are, respectively: /' , r' and g' for objects in the Local sample, and zgso , im and Veoe for those in the High-z sample. The ellipse 
marks the Petrosian radius (77=0. 2) in the B«f band and the solid horizontal lines indicate the equivalent to 5 kpc at the distance of 
the object. In the panels with the yU - r profiles, the vertical dashed line marks the position of the Petrosian radius. 
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Fig. 7. Left: Comparison our estimates of B -band luminosities (Mb ) of galaxies in the High-z sample, obtained with HyperZ, to 
those from COMBO-17. Right: Analogous to the left panel but for stellar masses (M* ). For both panels, the dashed lines represent 
a one-to-one relation, and the dotted, linear fits to the distributions of points. 



also applies to the vertical lines which mark the median values 
of the corresponding distributions. We see that there are some 
differences amongst the distributions of effective radii: for the 
Local sample the median value is Re/f =5.2+QA kpc, 4.3+0.2 
for the mid-z subsample (0.5< z <0.8 ), and 3.6+0.2 kpc for the 
far-z subsample (0.8< z <1.1 ). We recall here that no selection 
criteria were applied to the Local or High-z samples to put a 
limit on the maximum size of the galaxies, and so this differ- 
ence, though not very significant given the error bars, is not due 
to any obvious selection effect. According to this result, there 
is an increase by 44+20% in the R,,// (NUV) between z~l and 
z~0. If a Kolmogorov-Smirnov test is applied, the probability 
that the Local and mid-z samples of R^/f (NUV ) were extracted 
from the same distribution is rejected at 99.5% confidence level, 
while for the Local and far-z samples, this rejection is at a level 
>99.9%. If the same exercise is performed on images degraded 
to a common resolution (F =0.5 kpc, FWHMupc =1-5), the val- 
ues are not significantly altered. 

However, we do not claim that these figures are, on their 
own, proof of radial growth of the star-forming disc (or more 
precisely, growth in the distribution of A^f/V flux), because the 
galaxies in these samples do not share exactly the same ranges 
of luminosity and stellar mass, as seen in section |2] We thus 
present more reliable tests on this matter in next subsections. 

At this point it is convenient to explain why we select a max- 
imum distance, Dmax , of 60 Mpc in the selection of the Local 
sample, a limit that provides the largest possible sample while 
still having a negligible bias in size. To reach this conclusion 
we tested how the limit in distance affected the median value of 
Reff(BRf), Reff(B), in the Local sample. The median value of 
Reff(B) with Dmax =60 Mpc is 4.3+0.4 kpc (3 1 objects), whereas 
if that limit distance is 80 Mpc the sample is increased to 52 ob- 
jects, at the cost of an increase by ~35% in Reff(B) . To give a 
reference for considering the retrieved value of Reff(B) , we re- 
fer to Shen et al. ("2003), who estimated median sizes in a com- 
plete sample of galaxies, using different bands and estimators, in 
the local universe. The median r' -band luminosity in the Local 
sample is -20.95 mag. With this luminosity the disc-like objects 
in Shen et al. (I2003I I have Reff(r')—3.1 kpc. Shen et al. base 
their results on the extraction of radial profiles using circular 
apertures, i.e. they obtain circularized radii (r^). In contrast, our 
results are obtained taking as radii the semi-major axes of ellip- 
tical apertures (r^,) which fit the overall geometry of the objects. 
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Fig. 8. Effective radii of the rest-frame NUV flux distributions. 

As an approximation, = * where q is the axial ratio 
of the elliptical apertures. In the Local sample, the median value 
q =0.74 (corresponding to an inclination of 42deg), and so, the 
value by Shen et al. would correspond to ~4.3 kpc if they had 
used elliptical apertures. This result is in good agreement with 
ours, given that the effective radii measured in r' and g' are very 
similar (e.g. de Jong |1996] l. 

Another interesting test is to compare the spatial extension 
in the NUV to that measured in the B -band, which traces older 
stellar populations. In Fig.|9]we present the distributions of the 
ratio of the R^jj in NUVrf to that in B^f for galaxies in the red- 
shift bins considered. Also in this case there is some difference 
between the low and intermediate redshift samples. Galaxies in 
the High-z sample show median values of this ratio around unity 
(1.09+0.01 and 1.03+0.02 for the mid-z and fai--z subsamples 
respectively), while for those in the Local sample the median 
value is at 1.14+0.04. This means that local and higher redshift 
galaxies are slightly more extended in NUVrf than in Brf , but 
at z ~1 the difference is less. The null hypothesis with respect to 
the Local and far-z samples can be rejected at a level of confi- 
dence >99.9% (K-S test). 

A possible explanation of our findings could be the emer- 
gence of the bulge/pseudo-bulge in the galaxies as cosmic time 
evolves. In fact, all other structures in a galaxy being equal, a 
brighter bulge or pseudo-bulge would make the effective radius 
of the galaxy "shrink". This would be more evident in Bthan in 



12 



AzzoUini, Beckman & Trujillo: Radial Distribution of NUVFlux in disc Galaxies 



1 .0 - 



o 0.8 - 



0.6 - 



0.4 - 



0.2 - 



_z = 

..0.5<z<0.8 



_0.8<z<1 



O.OL^ 
0.0 



0.5 



R,„ NUV, 



1 .0 

..]/Re, 



1 .5 



2.0 



Fig. 9. Ratio of the effective radii in rest-frame NUV , NUVrf , 
and rest-frame B , B^f bands. 



NUV , as the populations in these structures are relatively old. 
We need to assume that there is not a parallel increase in the 
SFR in central regions that may cancel out this difference be- 
tween bands. Galaxies of earlier morphological types have, by 
definition, more prominent bulges. It could be argued then, that 
perhaps the apparent evolution in the ratio R^ff {NUV )IReff (B ) 
could be due to a difference in the distribution of morpholog- 
ical types, in which there would be a larger fraction of earlier 
types in the Local sample with regard to the High-z sample, and 
this could lead to a larger median ratio of effective radii between 
both bands in the former In consequence, we have tested the ra- 
tio R,ff (NUVrf )/Reff (Brf ) for local galaxies with 4< T <10 
(Sc to Sm), and the resulting median value is 1 .20+0.08, an even 
larger ratio. 

So, it seems that this increase in the ratio of effective radii 
is related to a progressive outwards migration of the SF rela- 
tive to the distribution of older stars in the discs. We have also 
compared the aforementioned ratios when we use equal reso- 
lution images (F =0.5 kpc, FWHMkpc =1-5), and the result 
stands, though the differences decrease a little, with median val- 
ues of the ratio R,ff{NUVRF)/Reff{BRF) of 1.13+0.03 (z~0), 
1.06+0.01 (0.5< z <0.8) and 1.02+0.01 (0.8< z <1.1 ). Again, a 
K-S test rejects the hypothesis that the Local and far-z samples 
come from the same parent distribution with high probability 
(confidence level >99.9%). In this way, it seems that the ratio 
Reff (NUVrf )/Reff (Brf) has slightly increased since z~l, by 
-10%. 



3.4. Median Intensity/Colour profiles 

In order to see with more clarity what the differences between 
the aforementioned effective radii in Brf and NUVrf mean, we 
also explore the median surface brightness profiles of the galax- 
ies in our sample. In Fig.[TO]we present median surface bright- 
ness profiles (// - r ) in NUVrf (blue squares) and Brf (red dots) 
for the galaxies in the Local sample (upper left) and at in- 
termediate redshifts (upper right, 0.5 < z <0.8and lower left, 
0.8< z <1.1 ). In the lower right panel are shown the NUVrf - 
Brf median colour profiles, for every redshift bin. In all panels, 
the radial coordinates have been scaled by the effective radius in 
Brf for each galaxy before taking the median. The dashed lines 
in the fi - r profiles give the dispersion (standard deviation). For 
the colour profiles, median values of the dispersion are 1.0 mag 
(z ~0), 0.8 mag (z -0.7) and 1.0 mag (z 



In Table [T] 7th column, we list the surface brightness level 
of the Icr fluctuations of the background, estimated in an an- 
nular aperture used with this aim in the production of inten- 
sity profiles and growth curves. These levels are fainter than the 
values of o-grcsec-^ for local objects, and brighter in the case of 
galaxies in the High-z sample, as the area of this annulus (and 
thus the number of pixels considered) is respectively larger, and 
smaller, than 1 arcsec^. The background subtraction is a signif- 
icant contributor to error in the profiles. We want therefore to 
estimate up to which radius the individual surface-brightness, 
and colour, profiles are reliable, taking into account this uncer- 
tainty in the background estimate. We are conservative in this 
regard, and following Pohlen & Trujillo (12 0061 see section 3.3 in 
their article), we adopt the following strategy to find these fidu- 
cial radii and surface-brightness levels. First, the corresponding 
Icfeflfi: fluctuations, converted to intensities, are added and sub- 
tracted to the median intensity profile in each band, giving the 
profiles ^{I + cThack, r) and yu(/ - cri,ack, f), respectively. Then, the 
fiducial radius is the largest for which the difference between 
these profiles, 6fi = |//(/ + o-back,r) - - o-back,r)\ < 0.2 
magnitudes. This is computed using the observed profiles, be- 
fore cosmological dimming, or any other corrections are ap- 
plied to the profiles (see below). In Fig. [10] we show the cor- 
responding fiducial radii and surface-brightness levels as ver- 
tical and horizontal dotted lines, in the same colour on the 
graph as the corresponding intensity profile. In Bs/rband, the 
fiducial radii/surface-brightness levels are at 3.1 Reff (Brf )/26. 1 
magjarcsec^ , and 2.7/25.9 and 2.8/25.4 for the Local, mid-z 
and far-z sample, respectively. In NUVrf the corresponding val- 
ues are 3.1/28.5, 2.4/26.6 and 2.5126.1 . So, the fiducial surface 
brightness levels are ~2.5 mag/arcsec^bhghteT than the IcTback 
fluctuations. 

For the colour profiles a similar approach is followed. 
Given the profiles of (5ju(r, band) in each band (NUVrf and 
Brf ), we obtain the fiducial radius as the largest for which 

(6^fi(r,NUV,-esi) + S^n(r,Bresi))^ < 0.2 mag. These radii are in- 
dicated as vertical dotted lines in the lower right panel of Fig. 
[TOl again in the same colour on the graph as the correspond- 
ing colour profile. The fiducial radii are at 3.1, 2.3 and 2.5 
Reff (Brf ), for the Local, mid-z and far-z samples. 

The profiles of the High-z sample shown in Fig. [TO] have 
been corrected for cosmological dimming, and a simplistic k- 
correction has been applied assuming a "flat" SED (fy - con- 
stant) and no inner dust absorption (Av=0 mag). These correc- 
tions are applied after the analysis of the reliability of the pro- 
files, described above, is performed. This is the reason why the 
lines of fiducial radius and surface-brightness do not cross at 
a point of the intensity profile for the mid-z and far-z subsam- 
ples, but at points 0.66 and 0.96 magnitudes fainter, respectively. 
Given these assumptions, the median value of surface brightness, 
ji, in NUVrf , at R -Reff , for galaxies at z~0 is roughly 1/5 of 
the value at z~l (A/i=1.8 mag). This is related to the significant 
decrease in the SFR the discs have suffered in the last ~8 Gyr. In 
the Bband this evolution is less significant, with A/i=1.0 mag; 
i.e. the surface brightness of discs (at R-Reff) in the B-band 
nowadays is ~40% of the value at z~ 1 . 

Fig. [To] helps us to understand Fig. [9] The Brf and 
A^t/Vfif profiles at intermediate redshift fall off exponentially, 
falling roughly "parallel" at all radii, because of the intense star- 
forming activity, which is largely dominating the appearance of 
the galaxies not only in NUVrf , but also in Brf . At z~0 the 
profiles also fall off, to a good approximation, with exponential 
rates, but the profile has a slightly steeper slope than the 
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Fig. 10. Median surface brightness (i-i-r) profiles (cosmological dimming corrected) of the galaxies in different redshift bins: 
z ~0 (upper left), 0.5< z <0.8 (upper right), and 0.8< z <1.1 (lower left), in B^f (red dots) and NUVrf (blue squares). The 
vertical dotted blue and red lines mark the radii up to which the NUVrf and Brf profiles are reliable, given the uncertainty 
in the estimate of the background, o-i,ack (i.e. the last radius with 6fi - {pil - (Tt,ack) - l^d + o-back)] < O.lmag, where cThack is the 
standard deviation of the background fluctuations in the typical annulus used to estimate it). The horizontal dotted lines are 
the corresponding surface-brightness levels at those fiducial radii. In the lower right panel are the corresponding median 
colour profiles (NUVrf -Brf )• There again, the vertical dotted lines in each color give the fiducial radii for the correponding 
color profiles {[d^fiB + S^I^nuvP < 0.2mag). In all cases the radial scale is relative to the corresponding effective radii, Re/f , 
of the objects in the Brf -band. In the yu - r profiles, the dashed fines give the standard deviation of the values of yu at each 
radius. The upward pointing triangles give the mean profile for galaxies with > 10 '"Mq , while the downward pointing 
triangles are the equivalent for galaxies with < 10 "'Mq . 



NUVrf in the inner parts of these median profiles. This differ- 
ence may contribute to make the ratio Rgjf (NUVrf )IReff {Brf ) 
at z~0 larger than the value at z~l. This central excess in 
brightness of the Brf profile at z ~0 could be the signature of 
older stellar populations being "piled up" in evolving bulges or 
pseudo-bulges. At z ~0, in NUV , there is also an excess in 
the central region over the disc profile. The central emission 
in Bsfband at z~0 is above the corresponding value at z ~0.7 
(A// ~ 0.6 mag), but only slightly above the corresponding value 
at z ~1 (Aju ~ 0.1 mag). In the NUVRF^and there is not a clear 
difference amongst redshift bins in this respect, though in the 
highest redshift bin (0.8< z <1.1 ) the central emission is the 
brightest. It must also be noted that the Brf and NUVrf profiles 
at z ~0.7 (upper right panel) show slight upbendings at small 
radii, similar to, but far less significant than those at z ~0, which 
we thus interpret as a hint of progressive evolution of pseudo- 
bulges in the explored range of redshifts. If we use images with 
equal resolution the shapes of the resulting median profiles do 
not differ significantly (the contribution of the bulges are more 
spread in radius, as expected), and also fit into the given inter- 
pretation. 

Also in Fig. \W\ we show the median profiles when galax- 
ies are segregated by their stellar masses: galaxies with > 



10 M0 are represented by upward pointing triangles, and those 
with M* < lO'^Moby downward pointing triangles. It can be 
seen that the shapes of the ju - r profiles when galaxies are segre- 
gated by mass are very similar to those of the whole populations, 
though there are minor differences may be seen. In Brf , more 
massive galaxies have somewhat higher levels of surface bright- 
ness than less massive. In NUVrf quite the opposite happens, 
median profiles of less massive galaxies, both in the Local and 
High-z samples, being somewhat brighter than those of more 
massive galaxies. This suggest a somewhat higher efficiency in 
SFR per unit area in less massive galaxies, a result that fits within 
the downsizing scenario. 

Another interesting feature may be found in the 
A^f/y profiles shown in Figure [10] Both for local and 
intermediate-redshift galaxies, the median profiles in this 
band are slightly "downbending", with the "elbow" of the 
profile located at R~Reff, something that is not seen in the 
B-band profiles. Here we recall that the intensity profiles of 
stellar discs may be classified, in general, in three types: Type 
I, in which the intensity decays exponentially with a certain 
scale radius, h (the slope of the fj. - r profile, m, is inversely 
proportional to this scale radius); type II, in which the profile 
declines exponentially with a certain h out to a break point, past 
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which a smaller value of h pertains, i.e. the descending slope of 
the profile becomes steeper; and type III, in which there is also 
an exponential decay in the intensity and again a change in the 
slope of the profile at a certain radius, but in this case the profile 
becomes shallower outside it (see, e.g. Erwin et al. .2008. for a 
review on related phenomenology). Stellar disc truncations are 
a subtype of type 11 profiles in which the break radius is located 
at what seems to be the "edge" of the stellar disc, beyond the 
extent of the spiral arms. 

The interpretation of the downbending feature in this figure 
is not straightforward, since the profiles shown here are the av- 
erage of profiles which have not been segregated according to 
type, and thus are mixed, supposedly in proportions similar to 
those found in field galaxies at z~Q and at z~l: ~60-50% of 
types II, ~30-40% of types I, and ~IQ% of types III (see e.g. 
Pohlen & Trujillo 2006 Azzollini et al. |2008bD . In any case, it 
seems that the A^f/V profiles show, commonly, a change in their 
slopes at some point in the discs, past which the drop in intensity 
becomes steeper 

Recently, it has been shown that truncated discs usually show 
radial colour profiles (u' -g' ) with a minimum value at the posi- 
tion of the break (i.e. it is a bluest point), and they become gradu- 
ally redder inwards and outwards from this point (Azzollini et al. 
I2008al Bakos et al. l2008l l. So, the colour profile has a character- 
istic "blue valley" shape. Even when the galaxies have not been 
segregated by the type of their intensity profiles (i.e. in types I, 
II, or III) the median colour profiles A^f/V -B presented in Fig. 
[TO](lower right) show a similar feature, most easily recognizable 
in the median colour profiles for galaxies at z ~0 (blue) and z ~1 
(red). This feature is clearly recognizable also when the FWHMs 
of the PSFs of the images are equal in both bands, as is the case 
when a uniform resolution (F =0.5 kpc, FWHMkpc =1.5) is sim- 
ulated. 

The local colour profile (blue squares) is clearly redder than 
those at z ~0.7 (green squares) and z ~1 (red triangles), by ~0.7 
mag, and the "valley" shape is more prominent, with a difference 
in colour AC~ 1 . 1 mag between the minimum and maximum (at 
the centre of the galaxies). In contrast, galaxies in the High-z 
sample show a more moderate minimum to maximum difference 
in colour of AC ~0.4 mag. We see that the minimum in colour 
takes place in the range \-\.5Reff ■ From these colour profiles 
we see how the central part is getting progressively redder in 
relation to the colour at the break. We interpret this as further 
evidence of secular bulge growth as cosmic time progresses. 

One question to be answered in the light of these results is 
whether the reported progressive rise of the inner parts of the 
median profiles, both in Brf and NUVrf bands, is due to a mere 
selection effect of the samples at different redshifts. As we said 
before, objects in the Local sample were selected to have mor- 
phological types in the range < T < 10, and those in the 
High-z sample, to have Sersic index n <2.5. We consider this 
as the best approach to have a sample of disc galaxies not bi- 
ased towards a particular Hubble morphological type, as already 
explained in Section |2l By definition, earlier type disc galax- 
ies have more prominent bulges, and so, if these were over- 
represented in the Local sample, with regard to the High-z sam- 
ple, this would cause a trend similar to that observed in the shape 
of the inner profiles with redshift. So, to settle this controversy 
it would be desirable to select the objects according to exactly 
the same criterion, for example the Sersic index, at all redshifts. 
Being short of that type of data, we have performed another test, 
to at least discard evidence for a bias within our limitations. We 
produced median intensity profiles for galaxies in the Local and 
High-z samples, dividing each of them into two more subsam- 



ples, termed as "earlier" and "later", according to their morpho- 
logical types. In the Local sample, the "earlier" subsample is 
composed of objects with < T < 4 (SOa to Sbc, 12 galaxies), 
and "later" of those within 4< T <10 (Sc to Sm, 19 galaxies). In 
the mid-z and far-z subsamples the "earlier" range corresponds 
to 1.25 < n < 2.5 (89 galaxies), and the "later" to < n < 1.25 
(149 galaxies). As would be expected, the median NUVrf and 
Brf profiles of the "earlier" subsamples show the trend to have 
brighter central parts with regard to the discs at lower redshift, as 
happens for the whole samples, yet amplified (i.e. with a larger 
contrast in brightness between the central part and the disc). 
Nonetheless the same effect, though somewhat less pronounced, 
is also present for the "later" subsamples. So, this effect does not 
seem to be due to a minor fraction of galaxies with prominent 
bulges in the Local sample, but is a common condition affect- 
ing, though in varying degree, most of the galaxies classified as 
disc-type, according to the prescribed criteria. 

There is also the important effect that variable dust extinction 
along radius may have on the shape of the Brf and particularly 
NUVrf , which so far has not been considered. We postpone the 
consideration of those effects until Sec. 13.91 where the effects of 
dust are treated with some detail, within the limitations of this 
work. 



3.5. Reff - B -Luminosity Relation 

A set of tests of particular interest for the study of disc growth is 
to compare the effective radii of these discs, as measured in dif- 
ferent wavelengths, to the luminosities and stellar masses of the 
galaxies. Here we focus on the R,,// as measured in the NUVrf , 
and so on the relation between the size of the star-forming disc 
and the stellar content of the galaxies, as a function of red- 
shift/time. 

In a first test, we study the relation between the Reff Va 
NUVrf and B -band luminosity {Mb ) of galaxies in different 
redshift bins (see Fig. [TT| for results). In the three panels we 
show Mh against R^jf in NUVrf ^ for galaxies in the Local sam- 
ple (z~0), mid-z subsample (0.5< z <0.8 ) and far-z subsample 
(0.8< z <1.1 ) respectively. The first thing to note is that there is 
significant dispersion in this relation at all redshifts. The vertical 
and horizontal dotted lines mark the median values of the dis- 
tributions for galaxies in the range -2Q<Mt, <-22 mag. The ob- 
jects have Z?R/r luminosities which gather around Mi, — 21 mag 
at all redshifts. The observed distributions have been fitted to a 
line with a least squares deviation fit, in an iterative "bootstrap" 
method, to get the most probable values for the slope and the y- 
intercept ordinate and their errors. These best-fit lines are shown 
as continuous in each panel, while the dashed lines in the second 
and third panels correspond to the best-fit line for the Local sam- 
ple relation, as reference. From visual inspection it seems clear 
that, at a given luminosity, Re/fis, smaller as redshift increases. 

In order to compare the radial extents of the SF of the galax- 
ies in a meaningful way, we can look at the values of the best- 
fit lines for the given relation, Rgff (NUVrf ) - Mt, , at a fixed 
value of luminosity, and plotted against redshift. This is pre- 
sented in Fig. [12] In it we show the best-fit R^jf value of the 
aforementioned relation at Mi,--2\ mag as a function of red- 
shift, in units of the corresponding value at z~0 (Local sample). 
The values taken from the best-fit lines are presented by empty 
diamonds, while the filled circles correspond to median values. 
The long gray error bars show the standard deviations cr, while 
the short black bars give the corresponding errors (cr/ V^, where 
A^ is the number of galaxies in the subsample). We see how the 
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Fig. 11. Effective radius of the rest-frame NUV flux distribution 
against the B -band luminosity (M;, ), for galaxies in different 
redshift bins: z=0 (Local sample, right panel), 0.5< z <0.8 (mid- 
z subsample, middle) and 0.8< z <1.1 (far-z subsample, right). 
The dotted vertical and horizontal lines mark the centroid of the 
distributions (median values), and the continuous line is the one 
that best fits the distributions. The dashed line in the middle and 
right panels gives the fit to the Local sample (left). 



median values and the best-fit values are in fair agreement. The 
best-fit points have been fitted to a line (continuous) which is 
forced to pass through the value at z=0. From this fit we obtain 
an increase by a factor 1.48+0.09 in Reff {NUVrf ) at Mt --21 
mag between z =1 and z =0. This analysis on the Br/t images 
yields a very similar growth factor (1.42±0.07). These numbers 
are not significantly altered when the images are degraded to 
the common resolution (F =0.5 kpc, FWHMkpc =1.5). If the 
median values are used, with the original resolution, the growth 
factor in NUVrf is somewhat lesser (1.35±0.07), and the same 
applies to Brf , with a factor 1.26+0.05, though these are still 
compatible with the best-fit values given the error bar. Again, 
very similar values are obtained if equal resolution images are 
used (F =0.5 kpc, FWHMkpc =1.5). For comparison, Trujillo 
& Aguerri (2004) found that galaxies with My < -18.5 show an 
increase in their V-band Ref/hy a factor 1.4+0.2 since z ~0.7. 
Barden et al. l2005l reported an increase in the same parameter by 
a factor 1.6+0.1 since z ~l (for galaxies with My <-20 mag). 

3.6. Reff - Stellar Mass Relation 

In the time interval between z~l and z~0, disc galaxies have 
undergone a significant, overall decrease in B -band luminosity 
(e.g. Rudnick et al. 120031 Barden et al. I2005I I. which is paral- 
lel to, and a consequence of the also reported decline in global 
SFR described in the Introduction. This situation adds some dif- 
ficulty to the interpretation of figures such as [TT] and [12] This is 
so because the observed evolution, with galaxies of a given lumi- 
nosity having larger R^ff 's at lower redshifts, has different inter- 
pretations depending on the relative contribution from changes 
in both parameters, R^ff and Mb . For example, in Fig.[TT]we can 
make the best-fit line for the far-z sample better match the Local 
sample either by moving it upwards, as would happen if galax- 
ies grew in R^ff but not in Mh , or shifting it leftwards, which 
would imply a pure decrease in M/, with time, but no change in 
the Reff of galaxies. 

To overcome the posed problem, we explore the relation be- 
tween Reff and stellar mass, M* . Working in terms of stellar 
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Fig. 12. Evolution of the relation Reff (NUVrf) - Mb. The 
empty diamonds give the value of ReffSt Mb=-2\ mag from 
the best fit lines in Fig.[TTl whereas the black dots are the corre- 
sponding median values of Reff . All measures have been scaled 
by the best fit value at z~0. The line that best fits the diamonds 
(continuous) gives an evolution in Rgff (NUVrf) ^.t the given 
luminosity by a factor 1.48+0.09 between z=l and z=0. The 
best fit line to the circles (median values of Reff , dashed line) 
represents an increase by 1.35+0.07 in the same redshift inter- 
val. The long-gray error bars give the standard-deviation (cr) in 
Reff (z)/Reff (z =0, while the short-black give the associated er- 
ror (cr/ -\/(N), where is the number of galaxies in the subsam- 
pie). 



mass is useful because it removes the evolution that is simply 
due to the ageing of the stellar populations. Here we recall that 
we want to test for evolution in the spatial distribution of SF 
in disc galaxies in relation to the process of stellar mass build- 
up, and so it makes sense to compare the radial extent of the 
NUVrf emission of galaxies with the same mass, at different 
redshifts. In Fig.[T3]we present the results of a test of this kind. 
There we can see the Reff (kpc) in NUVrf against the M* (Mq ), 
for galaxies in different redshift bins. Here we also see a signifi- 
cant dispersion in the relation, as in Fig.[TT] The dotted vertical 
and horizontal lines give median values in the represented pa- 
rameters for galaxies within 10^ <M* < 10" Mq . 

To the reported dispersion in the relations (Reff -Mb and 
Reff -Mi, ), in addition to their intrinsic origin, contribute the 
uncertainties in the distances to the objects, those in and 
Mb , and also, of some importance, the patchy distribution of 
NUV flux, dominated by conspicuous star forming regions. 

To test for evolution in the Reff -M* relation, here we choose 
to compare the best-fit values at = lO'" Mq , and the result 
can be seen in Fig. [14] As in Fig. [12] here we present the ra- 
tio Reff(z)IReff(z-Q) (best-fit values: empty diamonds) against 
redshift. From a linear fit (continuous line), we deduce an in- 
crease by a factor 1.18+0.06 in Reff (NUVrf) at = 10'° 

since z =l;i.e., a moderate change (<20%) in that quantity in 
the given period. If we measure the evolution in this relation by 
the median values of Rgff (filled dots, dashed line), the change is 
more significant (1.4±0.1), but this must be taken with caution, 
as the samples have slightly different median values of stellar 
mass, and Reff correlates with this parameter. 

Also, in Fig. [14] we show, as a dotted line, the evolution in 
the relation between the effective radius in Brf and stellar mass. 
The corresponding evolution is an increase by a factor 1 .10+0.05 
between z =1 and z =0, i.e. slightly less evolution than in 
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Fig. 13. Effective radius of the rest-frame NUV flux distribution 
against the stellar mass (M* ) for galaxies in the three redshift 
bins considered. This figure is analogous, in scheme and mean- 
ing of line styles, to Fig.fTTI 
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Fig. 15. Concentration of NUV flux for galaxies in diff'erent red- 
shift bins. C42 is defined as 5-log{R^()/R2o), where /Jgo and R20 are 
the radii which enclose 80% and 20% of the total flux of the ob- 
ject respectively. For the three redshift bins the median value of 
C42 is around 2.3-2.4 
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Fig. 14. Evolution of the relation Re/f (NUVrf ) - (see Fig. 
[T3] ). From the linear fits to this relation at M* = 10'" Mq (empty 
diamonds, continuous line), an increase inRef/hy 1.18±0. 07 be- 
tween z =1 and z =0 is deduced. If, instead, the median val- 
ues of Reff are used (filled dots, dashed line), the correspond- 
ing factor is 1.41 ±0.09. The dotted line shows the evolution 
for the Reff(BRf) - relation (1.10+0.05), from best-fit val- 
ues. The long-gray error bars give the standard-deviation (cr) in 
Reff {z)/Reff {z =0), while the short-black give the associated er- 
ror (cr/ yfiN), where is the number of galaxies in the subsam- 
ple. 

NUVhand. For comparison, Barden et al. (2005) reported prac- 
tically no evolution in R^ff in rest-frame V band for a fixed stel- 
lar mass between z ~l and z ~Q, in agreement with our re- 
sult within ~ 2cr. These figures also remain virtually unchanged 
when images of a same resolution (F =0.5 kpc, FWHMkpc -1.5) 
are used. 

3.7. Radial concentration and central emission of NUV -Flux 

So far we have presented results on the size of the galaxies as 
seen in NUVrf , and parameterized by the effective radii. In a 
step forward to characterize the SF distribution we have explored 
the radial concentration of the A^C/V emission in disc galaxies. 
This gives a measure of how relevant is the SF in the central re- 
gions of the galaxies in relation to that in the rest of the discs. In 



Fig.[T5]we show the distributions of the concentration parameter 
C42 in NUVrf band for samples in diff'erent redshift bins. The 
C42 parameter (Kent 1985) is defined as 5 ■ log{Rg,o/R2o), where 
R^o and R20 are the radii which enclose 80% and 20% of the to- 
tal flux of the object respectively. We find that the median values 
of C42 are around 2.3-2.4, and show no evolution with redshift 
within the error bars (Local: 2.35+0. 13, 0.5< z <0.8 : 2.28+0.05, 
0.8< z <1.1 : 2.39+0.05). The minor differences amongst red- 
shift bins remain minor when we analyse images at equal reso- 
lution, and are again not significant, given the error bars. 

But, what does this C42 parameter really mean? It is larger 
for those distributions of light which are spread over a larger 
area but at the same time have a larger fraction of the total flux 
accumulated near its centre. It may be better understood with 
some examples: for a disc with a flat surface brightness profile, 
C42 =1.51, for an exponential disc it is C42 = 2.75, and for a de 
Vaucouleurs profile (a Sersic with n=4), it is C42 =5.16 (Kent 
1985). So, the A^f/V profiles are between flat and exponential, 
being more similar to the latter, and their concentration does not 
seem to evolve in this period, at least by this parameter. 

The previous measure of A^t/V concentration depends only 
on the radial distribution of the emission in that band. But we 
also want to test how the SF is taking place in relation to the older 
stellar populations observed in the rest-frame B -band. In a first 
test, we have compared the fraction of the total NUV flux con- 
tained in an aperture with radius R -Rp/4, (Rpis the Petrosian 
radius, measured in the Bs/rband), for galaxies in the three red- 
shift bins under consideration, and the results can be seen in 
Fig- US] The median values of the fractions are 15.6+1.3% for 
the Local sample, 15.3+0.7% for the mid-z subsample (0.5< 
z <0.8) and 18.5+0.8% for the far-z subsample (0.8< z <1.1 ). If 
all of the NUV flux were contained in the Petrosian radius, and it 
were distributed uniformly in the area defined by that radius, the 
resulting fraction would be 6.3%, and so the results reflect a sig- 
nificant enhancement of the NUV emission in the central parts 
of the discs, which was expected from visual inspection of the 
images (also see Fig.fTOll. Also noteworthy is that the given mea- 
sure of central concentration of NUV finx remains virtually un- 
changed between z~0 and z~l, given the error bars. A K-S test 
shows that the null hypothesis can not be rejected with a confi- 
dence larger than ~98.5% in the best case, for the Local and far-z 
samples. The fractions obtained with images of uniform resolu- 
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tion (F =0.5 kpc, FWHMkpc =1.5) are not significantly different 
from the given. 

We have also performed this test with a smaller fraction of 
the Petrosian radius, Rp , to perform this test, Rp /IQ, and no 
significant difference amongst redshift bins is found either. For 
z~0 the fraction is 3.6±0.9%; for 0.5< z <0.8, 3.3+0.3%, and 
for 0.8< z <1.1 , 3.6+0.2%. At equal resolutions (F =0.5 kpc, 
FWHMkpc =1.5), the given fractions do not show any clear 
trend with redshift either: 3.1+0.5% at z~0 versus 2.6+0.1% for 
0.5< z <0.8 and 3.0+0.2 for 0.8< z <1.1 . 

The preceding measure of the fraction of A^f/Vflux con- 
tained in the "central" part of the galaxies is referred to the size 
of the galaxies in the Bband, through the Petrosian radius. The 
mean values of Rpait ~11 kpc for the Local sample, and ~8 
kpc for the High-z sample. So, given that these differences ex- 
ist, it is also interesting to define that "central" region without 
referring to any measure of the size of the objects whatsoever. 
To accomplish this, we also present results on the fraction of 
NUV flux contained within R<2.5 kpc, and this may be seen in 
Fig.[T7] Galaxies at z~l have a fraction of flux within R <2.5 kpc 
which is a factor <2 larger than that at z~0 (17.1+2.5% for z~0, 
24.1 + 1.8% for 0.5< z <0.8 and 31.1+2.2% for 0.8< z <1.1). 
Applying a K-S test, the null hypothesis for the local and z ~ 1 
samples can be discarded with a confidence >99.9%. These frac- 
tions are practically unaltered at equal resolution, and what we 
see here is that galaxies in the past had a larger fraction of 
A^t/yflux within a fixed radius {R <2.5 kpc), than nowadays. 
The difference is analogous, though reduced, if instead ofR <2.5 
kpc we use R < \ kpc, with a fraction of ~4% at z ~0 and ~6% at 
z ~1. But, given the error bars, the differences are not significant 
in that case. 

Another interesting question, and also related to the previ- 
ous one, is how intense is the Af/V emission in the central parts 
of the galaxies. We have chosen to give this measure in rela- 
tion to the brightness at the effective radius, Reff ■ In Fig. [18] 
we present the difference in surface brightness Aju, measured in 
rest-frame NUV , for the explored range of redshifts, and de- 
fined as Afi - fi{R = 0) - jj(R - Reff)- There we can see that 
this difference is more significant today than in the past, with 
median values of -2.0+0.2 (z ~0), -1.5±0.1 (0.5< z <0.8) and 
-1.2+0.1 mag/arcsec^ (Q.^< z <1.1). This is equivalent to an 
increase of brightness in the centre of the galaxy relative to the 
level at R^ff , from a factor of ~3 at z ~1 to ~6 at z ~0. In the 
B -band these differences are even more significant, with values 
of AyK of -3.1±0.2 (z ~0), -1.7+0.1 (0.5< z <0.8) and -1.5+0.1 



magjarcsec^ {Q.%< z <1.1 ), with an increase in relative bright- 
ness from a factor of ~4 at z ~ 1 to <20 at z ~0. For both bands, 
the hypothesis that the local and z ~1 samples are a realisation of 
the same parent distribution can be discarded with probabilities 
higher than 99.9%. 

If images of a same resolution (F =0.5 kpc, FWHMkpc - 1 -5) 
are used, the sense of these differences holds, though the values 
of A/i are somewhat lessened, and more similar to each other, 
as the central emission is smeared over a larger area when the 
images are degraded by convolution. This smearing effect is 
more prominent on images in the Local sample, as the differ- 
ences in spatial scale (F, FWHMkpc ) are larger within it. Under 
equal resolution, the values of A// in NUVrf still show evolu- 
tion {AjiiNUVrest) =-1.17+0.12 at z ~0, -1.11+0.05 at z -0.7 
and -1.05+0.05 magjarcsec^ at z ~1), although is less dramatic, 
and not significant enough under a K-S test (confidence level at 
81.4%). In contrast, in BpfXhty still show significant evolution, 
the samples being drawn from different distributions with a con- 
fidence higher than 99.9% {^^i{B) =-1.85+0.08 at z ~0 versus 
-1.21+0.03 at z~l). 

These results suggest that the slopes of the profiles, at least 
in the B band, have become steeper over time, or being more pre- 
cise, that the emission in the centres of the galaxies has greater 
contrast with respect to that in the discs. In fact, from the inspec- 
tion of Fig. [To] (left and central panels) it seems that the latter 
interpretation seems more adequate. There it can be seen that 
the slopes of the averaged profiles in B do not seem to vary sig- 
nificantly between z~l and z~0, but there is an abrupt increase 
in the emission towards the inner parts of the profiles for galax- 
ies at z~0, in contrast to the profiles at larger redshifts, which are 
flatter at smaller radii. 

The preceding results may seem, prima facie, and to some 
extent, contradictory. First we have that the concentration (C42 ) 
of the NUV emission does not seem to evolve in this range of 
redshifts. The same may be said about the fraction of flux in- 
side a certain fraction of the Petrosian radius, Rp (i.e. Rp/4-), for 
which if there is any evolution at all it is minor. We have also 
seen that the difference in brightness between the central parts 
of the galaxies and the disc (at R^ff ), is somewhat larger nowa- 
days than in the past. But, and here is the apparent contradiction, 
the fraction of flux inside a fixed radius (e.g. R <2.5 kpc) is larger 
at z ~1 than at z ~0 (this also applies, though not so significantly, 
if the central region is defined as that within R<1 kpc). The key 
to solve this "paradox" is that the first two measures of concen- 
tration are referred to the distribution of flux itself (in rest-frame 
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Fig. 18. Difference in rest-frame NUV surface brightness, Afi, 
between R =0 and R ^R^/f (Aju = n(R = 0) - n(R = Reff)), for 
different redshift bins. Median values of A// are -2.0+0.2 mag 
(z ~Q), -1.5+0.1 mag (0.5< z <0.8) and -1.2+0.1 mag (0.8< 
z<l.l). 



NUV or B ), and thus they are not affected by the fact that the 
objects are somewhat more compact (smaller) in the past, which 
the last measure shows as evident. Or as an alternative way to 
explain this, it seems that the centre of the galaxies has become 
slightly brighter in NUV in relation to the discs with time, but 
at higher redshift, the overall NUV emission, which had a some- 
what flatter profile, was also more compact (see also Fig. |8]and 
comments on this). 

Now it is useful to look again at Fig.fTO] to better understand 
these results. There we can see how the averaged A^C/V profiles 
(blue squares) for the Local sample show enhanced central emis- 
sion in comparison with the profiles for the High-z sample. The 
radial scale is relative to the effective radii, and so, in that fig- 
ure, the differences in compactness (i.e. size) of the objects are 
not evident. It is also important to note that the B profiles (red 
circles) also show enhanced central emission in the Local sam- 
ple. This seems to be related to the larger contribution of central 
bulges or pseudo-bulges to the profiles. So, it is doubtful that 
the slight increase in central NUV emission can be attributed to 
an enhancement of nuclear SF with respect to past times. Most 
probably this central NUV emission is due to intermediate age 
stars (with ages >100 Myr), at least in part. But a detailed analy- 
sis of the spectral energy distributions and stellar populations, 
in the central regions and in the discs, is necessary to solve 
this uncertainty. It must also be noted that the central excess of 
A^t/y brightness is roughly at the same surface brightness level 
at z~l and at z~0, and in Bthe central emission at z~0 is only 
slightly above the level at z~l. At the same time, the surface 
brightness of the discs has decreased significantly in both bands. 
So, it could be that the bulges/pseudo-bulges have "outshone" 
the discs over time, but, alternatively, it could also be that, as the 
discs have gradually faded, because of declining star-forming ac- 
tivity, the bulges have been left "exposed". At this point it is not 
possible to settle this dilemma. Nonetheless, we present some 
evidence in section 13.9.11 indicating that the second interpreta- 
tion may be closer to reality. 

3.8. Stellar Disc Truncations 

The galaxies in the High-z sample are taken from the survey on 
disc truncations presented in AzzoUini et al. ( I2008bl l. Thus, we 
have estimates of the "break" radii, Rb,- , for those galaxies com- 



mon to both samples which have truncated discs (137 galaxies). 
In the Local sample we find 13 truncations, and we have esti- 
mated their break radii. In both ranges of redshift, the band used 
to characterize the profiles is the rest-frame B -band. It is inter- 
esting to investigate the relation between the extension of the 
A^t/y emission (~SF), in relation to Rgr ■ This is because 7?^, is 
thought to be related to a threshold in the SF of the galaxies. We 
present results on this in Fig. \T9\ In it we can see the fraction 
of NUV emission contained in an aperture with radius equal to 
Rsr , for galaxies in the Local sample, and mid-z (0.5< z <0.8 ) 
and far-z (0.8 < z <1.1) subsamples. The distributions within 
the High-z sample (z ~1) are very similar, with median val- 
ues of 66+2% (0.5< z <0.8) and 63+2% (0.8< z <1.1). In the 
Local sample the fraction is significantly higher, with a median 
of 89±2%. The histograms of the mid-z and far-z subsamples 
reflect significant scatter, and it is interesting that the distribu- 
tions extend downwards to ~20%, implying that there are discs 
in which most part of the NUV flux outside the break radius, 
contrary to what happens for the histogram at z ~Q. Here, again, 
using images with equal resolution does not have significant con- 
sequences for the results. 

So, it seems that local truncations include most of the 
NUV flux, and hence, presumably also the star-forming activ- 
ity. This fits into a scheme for the formation of truncations in 
which the suppression of SF outside the break plays a major role. 
Nonetheless, the fraction of A^t/y flux within 7?^,. at z ~1, though 
high, seems to challenge this picture, if we take the emission in 
that band as a tracer of on-going SF. But several things must be 
said before reaching any conclusion. First, it is known that the 
profile classifications are somewhat less reliable at higher redsh- 
fit due to poorer signal-to-noise ratios. For example, in ATB08b, 
and from simulations of artificial galaxies, it was estimated that 
~30% Type 1 galaxies of magnitude zaso ~23 are wrongly clas- 
sified as type 11. Second, the NUV flux is not a perfect tracer 
of SF, as has been said, and stars with ages of a few xlO^ Myr 
also contribute to the emission at those wavelengths. So, stars 
recently formed inside the break which have migrated outwards, 
as in the models by Roskaret al. ( 2008all20 08b). could also con- 
tribute to the A^[/y emission outside the break. Nonetheless, it 
could also be possible that there might be somewhat more SF 
outside the breaks at higher redshift than in the local universe 
perhaps as a consequence of a more important merger activity in 
past times, and this could explain the reported difference, or part 
of it. 



3.9. The Attenuating Effect of Dust in the B and NUV bands 

In the preceding sections we have presented results on the radial 
distribution of NUV flux, analysing several parameters that char- 
acterize these distributions and also median surface brightness 
and colour profiles. We are ultimately interested in the spatial 
distribution of the SF, but as we have said, the effect of dust may 
play a significant role in the shaping of the A^f/y profiles, and so 
hinders an interpretation of the given results in terms of SF. Here 
we show some tests we have performed in order to estimate the 
effects of dust attenuation on the results. 

In Boissier et al. (I2004I I. an analysis of the radial FIR to FUV 
ratio in a sample of 6 local disc galaxies was used to provide a 
recipe for dust attenuation in the UV (at 2000 A) with radius, as 
a function of the B -band luminosity of the galaxy, and its incli- 
nation (see equations 5, 6 and 10 in their paper). Assuming the 
dependence of attenuation on wavelength (k^) given by Calzetti 
et al. ( I2OOOI 1, we can get an estimate of the extinction as a func- 
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Fig. 19. Fraction of NUV^ux within the "break" radius (mea- 
sured in B/tf ) for galaxies in the Local and High-z samples. The 
classification of profiles and estimates of Rsr for objects in the 
High-z sample is taken from Azzollini et al. 12008b] 



tion of radius in the bands explored (NUVrf and Brf ), Abandir). 
The resulting extinction profile significantly obscures the cen- 
tral parts of the galaxies {Amuv{Q) ~2.5 mag), while it becomes 
a minor effect beyond R ~2Reff {A^uviR > '^Reff) <0.4 mag). 

In Fig.|20]we show the median surface-brightness and colour 
profiles of Fig. (TO] after allowing for the dust extinction in the 
NUVrf and Br/t bands. We have used median values of inclina- 
tion and absolute B -band magnitude for each sample to compute 
the extinction profiles and these have been subtracted from the 
stacked profiles. We can appreciate how the yu - r profiles be- 
come steeper in the inner parts after this correction. The previ- 
ously noted hints of profile downbending in the A^J/Vr/t profiles 
are washed out in the Local sample, and decreased significantly 
for the higher redshift samples. The colour profile for the Local 
sample still has a valley-like appearance, though for the mid-z 
and far-z sample this shape has disappeared, and instead we get 
rather flat colour profiles (out to the fiducial radius), with a pos- 
itive slope (i.e. the median profiles become redder to the outer 
parts of the galaxies). The colour profile of the Local sample is 
now not so red in the central parts, as it was before the correc- 
tion. 

Can we trust these extinction corrected profiles? Before in- 
terpreting these figures, it must be emphasized, though, that the 
extinction curve by Boisser et al. was derived from the analy- 
sis of a very small sample of local galaxies. So the reliability 
of the applied correction is even less at intermediate redshifts 
than it is in for the Local sample, as there might be evolution- 
ary changes in the radial distribution of dust attenuation in this 
time interval. This is probably the case as as the SF distribution 
is quite different at intermediate redsdhifts than currently, and 
consequently the dust distribution is also likely to be different. 
Nonetheless, this result indicates that in the progressive redden- 
ing of the inner parts of the profile since z ~ 1 , the ageing of the 
stellar populations must play a major role. 

Also in Fig.|20l we can see how the rise of the central parts of 
the ^ - r profiles at lower redshift, in comparison to the profiles 
at higher redshift bins, still holds. We have not modelled any 
possible evolution in the radial distribution of dust extinction, 
and so this difference was not expected to be altered after this 
correction. Only a detailed analysis of the dust extinction profile 
in each galaxy, from the modelling of the SEDs of the galaxies, 
could provide a more accurate picture. 



But we have also explored the net effect of dust on the effec- 
tive radius and other parameters of the radial profiles commented 
in preceding sections. With this aim, we have taken the median 
profiles in Brf and NUVrf for the Local sample, subtracted the 
corresponding radial curves of extinction (Boissier et al. I2004I I. 
and measured the parameters in both sets of profiles, those cor- 
rected and not corrected from extinction, for comparison. 

The effective radii seem to be significantly affected by dust 
attenuation, and after dust correction, these radii decrease, in 
Brf and NUVrf , to 61% and 64% of the non-corrected values, 
respectively. This is a direct consequence of the assumed radial 
extinction curve, which implies higher absorption in the central 
parts of the galaxies. However, the differential effect that dust 
has on ^e//in both bands could only account for part (<6%) of 
the measured difference between the effective radii, by which the 
in NUVrf is ~10% larger than in Brf ■ 

The concentrations (C42 ) are also affected by the dust cor- 
rection. The dust corrected value in NUVrf , is C42 =2.52, ver- 
sus 2.25 without correction. The new value is closer to the value 
for an exponential profile, but still below it (C42 =2.75). Perhaps 
more interesting are the results with regard to the fraction of 
flux contained in a given fraction of the Petrosian radius or in a 
fixed radius. On one hand, the former is practically unaffected, 
the fraction of NUV^ux contained within Rp/4 being 19.7% 
and 21.6% before and after the dust correction. This could be 
anticipated, as the extinction curves are very similar in shape be- 
tween these two bands, as made evident by the little differential 
effect they have on the effective radii, commented above. On the 
other hand, the fraction of NUVfiux within a fixed radius (2.5 
kpc) does show a significant difference, being 20.6% before the 
correction, and 41.1% after it. This is so because the profile is 
particularly attenuated in the inner parts of the galaxies. 

The difference in surface brightness between the centre and 
Reff^s also affected by the dust correction. Without dust cor- 
rection, Afi{NUV) =-1.5 magjarcsec^ , and after correction it 
becomes -1.7 mag/arcsec^ (increasing). In contrast, in Bband, 
it goes from -2.4 to -2.2 mag/arcsec^ (decreasing). This differ- 
ence in the sense of the variation is also due to a larger dust 
absorption in NUV than in B (the reader must bear in mind that 
the definition of A/i(band) contains R,,// (band), which is also 
affected by the dust correction). 

We have gone one step further, and performed the same anal- 
ysis described here, this is, measuring the parameters studied 
before and after dust correction, but using the median NUV and 
B profiles at z ~1 (far-z sample). Combining these results with 
those obtained using the z ~0 profiles, it is possible to estimate 
the contribution that dust absorption makes to the observed vari- 
ation of these parameters in this time interval. These net effects 
are: 1) Reff va B increases by 17% since z ~1 due solely to the 
effect of dust. This value is similar to the evolution found (an 
increase by 10% at a fixed mass), and so perhaps the reported 
evolution could be mostly due to dust. 2) In contrast, dust does 
not seem to have a net effect on the evolution of Re/ j- {NUV), 
and thus the reported evolution (an increase by ~20%), seems to 
reflect a change in the distribution of young stars in the galax- 
ies. 3) As a consequence of 1) and 2), without the effect of 
dust, the ratio Reff{NUVRF)IReff {Brf) would have increased 
by ~25% since z ~1, above the reported increase by 10%. 4) 
The Concentration, C42 , in A^t/VR/r seems to decrease by 0.21 
points since z ~ 1 . Taking this into account, C42 would have in- 
creased slightly since z ~1, were it not for the dust effect (the 
direct measure did not show a clear evolution of this param- 
eter). 5) The dust makes the fraction of flux inside the cen- 
tral 2.5 kpc decrease by 12%. The reported evolution is a de- 
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crease by ~55%, and so the evolution of the fraction of SF 
would be somewhat lower, closer to 40%. 6) The fraction of 
flux inside Rp/4- increases by 17% due to dust, and so, as we 
measure no evolution in this parameter, this suggests that the 
fraction of SF inside that radius has actually decreased, since 
z ~1, by that relative amount. 7) The net effect of dust on the 
difference in brightness between centre and R^ff is not big in 
either band. The variations attributable to dust would be -0.08 
mag /arcsec^ and +0.12 mag/arcsec^ in B«f and NUVrf since 
z ~1. Accounting for these trends, the real evolution of A/u in 
Buf and NUVrf would be, respectively, slightly less and slightly 
more significant than those found. 

After these tests, we conclude that the effect of dust absorp- 
tion on the profiles is significant, in absolute terms, for most of 
the studied parameters, especially for the effective radii, the ex- 
ception being the fraction of A^C/V flux inside a fraction of Rp . 
Nonetheless, being short of a prescription for the evolution in the 
radial distribution of dust, we cannot estimate in a reliable way 
the effect of the dust absorption on the evolution of the explored 
parameters (e.g. R^fj (NUV)). Assuming that the extinction law 
for the local galaxies (from Boissier et al. 2004) also applies at 
z ~l, which is at the very least dubious, it seems that the ob- 
served increase of R,,// {NUV )/Reff (Brf ) with time is a lower 
bound to the underlying evolution, once dust absorption is ac- 
counted for. On the other hand, the progressive decrease in the 
fraction of SF within the central 2.5 kpc would not be so large 
as direct measurement suggests, though still significant. Finally, 
there are two evolutionary trends whose reliability is reinforced 
after these tests. First, the central parts of the galaxies are red- 
denning largely due to ageing of the stellar populations (becom- 
ing brighter in Bwith respect to NUV). And second, the effec- 
tive radius of the SF seems to be increasing with time. Instead 
of the last assertion, perhaps the reader would have expected the 
more straightforward "the SF is migrating outwards". But before 
stating that this is what is actually going on, we should read the 
next subsection. 

3.9.1. SFR Radial Profiles and their Integral in Time 

Based on the Bpf and NUVrf profiles presented in Figs. [TOl and 
l20l we have also explored the stellar mass density profiles of 
our samples of galaxies. In particular, here we compare the esti- 
mated present profile of surface mass density with the accumu- 
lated mass profile, obtained from the time integral of the SFR 
radial profiles since z ~1. We rely on the chosen dust correc- 
tion (Boissier et al.'2004!) to estimate the SFR profiles from the 
NUVrf surface brightness profiles, but it must be taken into ac- 
count that this correction was derived from a reduced sample 
of local objects and we apply it to objects expanding in a wide 
range of redshifts. Furthermore, we sum the SFR over time to 
produce a profile of accumulated stellar mass, and compare it 
with the current stellar mass profile. But the NUVrf profiles at 
intermediate redshifts used do not exactly correspond to the sus- 
pected progenitors of the galaxies in the Local sample. We re- 
call that the Local sample has a median value of stellar mass 
log(M0)=lO.l, while in the High-z sample it is log(M0)=9.6. 
This difference is above the reported evolution in stellar mass be- 
tween z ~l and z ~0 (<50%). So we warn against basing strong 
conclusions on this test. Nonetheless, the results are interesting 
enough to deserve mention. 

The local Br;? profile provides the best proxy, available to 
us, for the radial distribution of mass, and we have converted it 
to a surface mass density profile (S'^*(r), in units of Mqpc^^), 
using Eq. (1) in Bakos et al. ( 120081 ). The mass to luminosity 



ratio in g' band, (M/L)g>, is estimated using results in Bell et al. 
( 120031 see Table 7), and assuming a Salpeter (1955) IMF. We 
use median rest-frame g' and r' profiles, without correcting for 
dust, to estimate the surface density mass profile. We proceed in 
this way because the recipe for (M/L)g> given in Bell et al. does 
not account for such extinction corrections. 

With regard to the profiles of surface SFR density, 
i.e. S'^^^(r), given in units of MQyr^^ pc^^, we use the 
NUVrf profiles within each redshift bin, after dust correction. 
These are transformed to the sought-after profiles via Eq. (1) in 
Kennicutt ( 119981 ) (the Salpeter IMF is also assumed for this ex- 
pression). These surface mass density and surface SFR density 
profiles are re-scaled into kpc, by means of the median values 
of Reff {Brf ) (in kpc) within each redshift bin (see Section 3.3). 
Then the surface SFR density profiles are linearly interpolated 
and integrated in time to produce profiles of the stellar mass pro- 
duced at each radius (in kpc) since z ~l (accumulated mass). 

In Fig. |2T1 upper panel, we show the radial profiles of SFR 
surface density for the three redshift samples. We can see that 
the surface density of SFR has decreased since z ~l, and that 
these decrease has been somewhat greater in the inner parts of 
the galaxies (within the inner ~2.5 kpc). And it is remarkable 
that the profiles are roughly parallel beyond R ~4 kpc. These 
trends would cause the effective radii of NUV to increase with 
time, as we had reported. The median effective radii, in kpc, of 
the dust-corrected NUVrf emission can be seen as vertical lines 
in the same colours as the corresponding curves of SFR density. 
Also note that the innermost part of the local profile has an up- 
ward bend which is probably due to the contribution from not so 
young stars (with ages around ~ 10^ yr), and not an indication of 
recent SF. This central "peak" in NUVrf (from which the SFR 
profile is derived) is responsible for the increase in the value of 
Afj. since z ~l that we have also reported. 

The blue dotted curve in the central panel gives the amount 
of stellar mass (density) produced at each radius since z ~l- 
Also in the same panel, and shown as a continuous red curve 
is the present stellar mass density, derived from the stacked 
Brf profile (z ~Q), as explained above. If the accumulated stellar 
mass is subtracted from this profile we obtain an approximation 
to the mass profile at z ~1, shown as the dashed red line. The 
dashed and continuous vertical lines show the effective radii of 
the Br/t profiles at z ~Q and z ~l- The backwards extrapolated 
mass profile also shows hints of a pseudo-bulge, as the local one 
does, and so this suggests two scenarios. First, it could be that 
this feature was already largely in place at z ~1, but was out- 
shone by the emission from young stars. This emission from 
young stars contributes significantly to the emission in Brf at 
z ~l, and so, this pseudo-bulge structure would not be recog- 
nizable in the z ~1 profile in this band. This, however, would 
need confirmation from high-resolution observations in redder 
wavelengths at that redshift. Alternatively, it could also be that 
this pseudo-bulge structure we see in the local Brf profile were 
the by-product of some mechanism of stellar mass redistribution 
within the galaxies in this time interval. 

It is also noteworthy that the maximum variation in stellar 
mass is not at the centre. This is easier to appreciate in the lower 
panel, where we show the relative increase in mass since z ~1, 
according to this test. The mean increase between R= Q and 
7? =8 kpc is ~120%. But the maximum in the relative increase 
reaches ~300% at /?~ 1.5-2 kpc, possibly indicating a progres- 
sive dearth of the SF activity in the central part of the galaxies. It 
could be that a higher efficiency of the SF in the central parts 
(perhaps due to a higher gas density), had exhausted the gas 
reservoir more quickly in that area. This would be at a greater 
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rate than that at larger radii, where a decline of SF with time 
also seems to happen, though not so marked. This indicates that 
either the parameters that regulate the efficiency of SF vary with 
radius, or that the gas replenishment is more effective in the outer 
parts of galaxies. An alternative, though admittedly less proba- 
ble explanation, is the nuclear-activity-induced quenching of the 
SF. 

Interesting as this may be, we need to present the caveats 
affecting this result. First and foremost, it is based on rather 
inaccurate dust extinction corrections which neglect possible 
evolution in the radial distribution of dust. Second, the local 
Bfif profile is not the most reliable proxy for the stellar mass 
density. Similarly, the NUVrf profiles are also "contaminated" 
to some extent by the contribution from not- so-young stars. 
We must also note that the stacked profiles are for collections 
of galaxies which have similar stellar masses, luminosities and 
morphological types, but within a range, and they are not se- 
lected to represent the possible evolutionary path of the popula- 
tion of disc galaxies (i.e. trace their mass evolution). And last, 
but not least, we have ignored the effect of stellar mass redistri- 
bution within galaxies, and amongst galaxies (through mergers). 
These mechanisms most probably would make the true surface 
mass density profile at z ~ 1 differ from the given approximation, 
which is based solely on the SFR curves. So, the main interest of 
this exercise is to test whether the SF activity in the inner parts of 
the discs could be, alone, responsible for the pseudo-bulge-like 
structure seen in the local B«f median profile. And the answer 
seems to be "no". This test suggests that some other mechanisms 
must be collaborating in this phenomenon. 

We have also checked that the profiles shown in Fig. 
l2n are basically unaffected when equal resolution images are 



used. Besides, the fiducial radii in Buf and NUVrf are be- 
yond the shown radial ranges, and the limit is imposed by the 
NUVrf profile at z ~ 1 (marked as a dash-dotted black vertical 
line in the panels in Fig.l^TTl. 

4. Discussion 



Barden et al. ( I2005I I found that the B -band effective radii of stel- 
lar discs, at a given stellar mass (M* > lO'^M© ), do not evolve 
between z~l and z~0 (although this result seems to be mass 
dependent, see Trujillo et al. 2007). Given that galaxies have in- 
creased their stellar mass content in this period (by ~3Q% ac- 
cording to Rudnick et al. I2003I I. these findings imply that indi- 
vidual disc galaxies have grown in size {Re/f ), so as to keep their 
surface mass densities roughly constant. The absence of evolu- 
tion in stellar mass density is at odds with some previous theo- 
retical expectations. Mo et al. 1 199 81 proposed a model in which 
the size of a disc of fixed mass is to scale with time in proportion 
to the virial radius of the dark matter halo. In this model, this is 
the prescription for the evolution in effective radius, Re/f , at a 
fixed mass: 



Reffiz) ^ Reff(O) 



H{z) 



H(Q) 



-2/3 



(1) 



where Ref/iO) is the effective radius at z -0 and H(z) is 
the Hubble parameter as a function of redshift. According to 
this, stellar discs at z ~ 1 would be a factor 2 denser than they 
are nowadays, in contradiction with the cited observational evi- 
dence. 

When the expected disc size-stellar mass distribution is com- 
puted, including the evolution in the internal structure of dark 
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Fig. 21. Upper: median SFR surface density profiles at differ- 
ent redshifts, derived from the dust-corrected NUVrf profiles. 
The vertical lines are at the the corresponding median values 
of Reff (NUVrf ), corrected from the effect of dust. In all pan- 
els, the vertical black dotted line marks a 2.5 kpc radius, and 
the vertical dot-dashed line marks the fiducial radius for the 
A^f/VRf profile at z ~1, which sets the limit for the reliability 
of the test, because of errors in the estimate of the background. 
Central: median stellar mass surface density profile at z ~0, de- 
rived from the dust-corrected Br/t profile in Fig. |20] (continuous 
red curve). The blue dotted curve is the integral of the stellar 
mass born at each radius since z ~ 1 , from the linear interpolation 
and time integral of the SFR curves in the upper panel. The red 
dashed line is the result of subtracting this "accumulated" mass 
profile from the present stellar mass profile, i.e. a backwards ex- 
trapolation of the present day mass profile to approximate that 
at z ~ 1 . Lower: Relative increase in stellar mass since z ~ 1 as 
a function of radius, from the ratio of the current stellar mass 
density profile to the backwards extrapolated profile at z ~ 1 . See 
the details in the text (Section D.9.1t . 



matter halos and the adiabatic contraction of the dark matter by 
the self-gravity of the collapsing baryons, the predicted evolu- 
tion in the mean size at fixed stellar mass since z ~1 is about 
15-20% (Somerville et al. I2008I I. in reasonable agreement with 
the aforementioned observational constraints. 



Our results on the evolution of the effective radii of galax- 
ies in the B -band agree fairly well with the results by Barden 
et al. 2005, and so, they are also consistent with the model by 
Somerville et al. 120081 But the effective radius in NUV seems to 
grow somewhat more than in B -band in this time interval (rel- 
ative to Brf , by ~10%). We have seen that at least part of this 
difference is due to a change in the shape of the profiles in both 
bands. It seems that the profiles have changed from "flat" in the 
central parts of the galaxies at z~l, to more "cuspy" at z~0, an 
evolution which is stronger in B -band. We attribute this change 
in shape to the interplay of two processes: a) a progressive pil- 
ing up of stellar mass in the central parts, i.e. the signature of 
growing bulges and pseudo-bulges in discs, and b) to the age- 
ing of stellar discs, due to a decreased star-forming activity per 
unit area, leading to a progressive dimming of the discs in the 
NUV and B bands. From our results the relative contributions of 
these mechanisms to the final observed evolution of the shape of 
the profiles cannot really be disentangled. The bulges/pseudo- 
bulges, which are easily seen in the z~0 profiles, could already 
have been in place at z~l, and not have evolved significantly 
in this interval, but could be outshone by bright, actively star- 
forming discs at z~ 1 . The fact that central brightness in the mean 
NUV and B profiles has not evolved significantly over this range 
of redshifts, as seen in Fig. \W\ while the brightness of the discs 
has significantly decreased, argues in favor of this interpretation. 
Furthermore, when we try to reproduce the mass density profile 
at z ~1 from the present curve (^Br/t profile), and the time inte- 
gral of the SFR profiles { ~NUVrf profiles), though applying less 
than fiducial dust corrections, leads to a past mass profile where 
this pseudo-bulge feature is already present. On the other hand, 
the strong evolution in colour in the central parts {R ~0), being 
much redder nowadays (even after dust-corrections), favours the 
interpretation that there might have been substantial secular evo- 
lution of bulges/pseudo bulges in the last 8 Gyr 

The concentration (C42 ) of the NUV profiles does not seem 
to vary significantly in this time interval (0< z < 1), and the val- 
ues of this parameter indicate that globally the profiles are close 
to exponential, but somewhat "flatter". The fraction of flux en- 
closed within a fraction of the Petrosian radius (measured in B ) 
seems also not to evolve. But the fraction of flux within a fixed 
radius, given in kpc (e.g. R <2.5 kpc) shows significant evolu- 
tion, galaxies at z~ 1 having a larger fraction of NUV flux within 
that radius than their peers at z~0. This could indicate that the 
star-forming area in them was more compact in the past, but the 
distribution in relation to that at longer wavelengths (B ) does not 
change greatly. An evolution of this kind suggests that the stellar 
discs are growing from the inside outwards. 

The difference in central brightness between the centre and 
7?(,//has also evolved, the objects at z~0 having their central 
parts brighter in NUVwd especially in B by a larger factor than 
the discs, in comparison to galaxies at z~l. This trend in the 
B band seems to be largely due to current bulges/pseudo-bulges 
being more prominent with respect to the discs than in the past. It 
remains to be answered whether the similar evolution observed 
in NUV is due to a parallel boost of nuclear star formation in re- 
cent times, or to the contribution of numerous middle aged stars 
in those central areas, or to some intermediate situation. The me- 
dian colour profiles for galaxies at z ~0 (Fig. [TOl ) show signif- 
icant reddening at small values of R , but this, of course, does 
not suffice to state that all A^f/V emission in the central parts 
of galaxies is due to intermediate age populations. A detailed 
study of the composition of stellar populations, via the analysis 
of SEDs, in the central areas of discs would be required to clarify 
this issue. 
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In this work we have also explored the relation between 
the radial distribution of NUV^nx (a proxy for SF), and the 
phenomenon of stellar disc truncations. The mechanism behind 
these truncations is not fully understood yet, though several 
models have been proposed. Most of them can be grouped in two 
groups: a) models related to angular momentum conservation in 
the proto-galactic cloud (van der Kruit 1987 ) or angular momen- 
tum cut-off in cooling gas (van den Bosch 2001 ); and b) models 
which appeal to thresholds in star formation (Kennicutt 1 19891 1 
or star-forming properties of distinct ISM phases (Elmegreen 
& Parravano 1994; Schaye 2004; Elmegreen & Hunter 2006). 
More recently, N-body simulations in which both star formation 
processes and secular dynamics are included, have been devised 
in an attempt to better fit observational results (e.g. Debattista 
et al. 2006). In a further elaboration of these ideas, in Roskaret 
al. (|2008a| ) a model was presented in which the breaks are the 
result of the interplay between a radial star formation cut-off 
and a redistribution of stellar mass by secular processes (see also 
Roskaret al. |2008b] l. In this model, stars are created inside the 
break (threshold) radius and then move outwards due to angular 
momentum exchange in spiral arms. Other recent models ex- 
plain the occurrence of breaks in the mass density profiles as the 
result of other secular processes. Foyle et al. (2008) presented 
a model in which angular momentum redistribution inside the 
break was the cause for it to form. In a different approach, in 
Bournaud et al. (2007) it was shown that the disruption of large 
stellar "clumps" (with stellar masses in the range 10** - 10^ Mq ) 
in primordial discs (observed as "clump cluster" galaxies at high 
redshift) could give way to single or double exponential discs 
(i.e. discs with a break in their profile). 

We have presented results on the fraction of NUV^ux con- 
tained within the break radius, Rgr , for galaxies in the Local and 
High-z samples. This fraction is relatively high, but there is sig- 
nificant difference between the local and intermediate-redshift 
results, the median value at z ~0 being ~90%, while at z~l it 
is >60% (in some cases it is as low as ~20%). The values at 
intermediate redshifts seem perhaps too low, assuming that this 
NUV emission traces the star-forming activity, to fit into mod- 
els in which the suppression of star formation outside the break 
radius is part of the mechanism behind the formation of a trun- 
cated disc. Being strict, one would expect a distribution simi- 
lar to that of the Local sample. But other factors must be taken 
into account before extracting hasty conclusions from this re- 
sult. First, the NUV is not a perfect tracer of recent SF, and so, 
even if the break could mark the extension of the area where 
SF takes place, stars outside this break, of intermediate age, 
could also contribute significantly to the total emission of the 
galaxy in that band. For example, in the akeady cited model by 
Roskar et al. ( 2008a,, 2008b ), it is proposed that a mechanism 
of radial migration of stars due to secular processes may drive 
stars away from the zones where they form, supposedly inside 
the "break". Furthermore, from visual inspection we have seen 
that in some cases there are few, but powerfully emitting regions 
outside Rb, which account for reduced fractions of NUV flux en- 
closed within R<Rbi ■ It could also be possible that there was 
some contribution to the A^t/V emission outside the break due 
to mergers. In any case, this result does not suffice to decide 
what is the mechanism behind disc truncations, whether models 
which rely on angular momentum conservation / cut-off in the 
proto-galactic cloud, or those which postulate that the break is 
the signature of a threshold for SF in the disc. At most, it can be 
said that the fraction of SF outside the break at z~ 1 seems to be 
rather high, even given the limitations of NUV flux as a tracer 
for the star-forming activity. 



We have seen that the median B and A^t/V profiles have de- 
creased their brightness at R -R^ffhy ~6Q% and ~80% respec- 
tively between z~l and z~0. Although these variations have 
been deduced assuming simplistic k-corrections, this result is 
consistent with a significant decline in the SFR per unit area in 
stellar discs in the last ~8 Gyr. The decline in NUV , and in 
B , has been similar for galaxies within different ranges of mass 
(M* < 10'°Mq and > lO'^M© ). 

Summarizing, stellar discs, as seen in the NUV , seem to 
have undergone significant evolution on several fronts, and this 
could be succinctly described as follows: a) the surface bright- 
ness of the discs has decreased to roughly a 1/5 of the median 
values at z~l; b) NUV discs have grown, at a fixed stellar mass, 
only moderately (by <20%); c) the ratio of the effective radius 
in NUV to that in B has increased by -10%; d) the NUV profiles 
are rather flat in the inner regions, showing signs of "trunca- 
tion" in many cases, as made evident in averaged profiles for 
the whole samples at different redshifts; e) a significant excess 
in emission (NUV and B) at the centre of the discs has devel- 
oped over time (between z~l and z~0). What does this informa- 
tion tell us about the distribution of SF? If the relation between 
A^f/y emission and SFR were strictly proportional, it could be 
said: a) the star-forming activity decays with radius according to 
a double exponential law, with a larger scale radius (i.e. a shal- 
lower slope of the profile) in the inner parts of the discs; b) the 
SFR per unit area in discs has decreased significantly in the red- 
shift interval 0< z <1; in this same period c) the bulk of this 
star-forming activity has shifted to larger radii within the discs; 
and d) the SFR in the bulges/nuclii has grown in relevance over 
that in the discs. Nonetheless, NUV flux is not a perfect tracer of 
this activity, even when it provides a fair approximation. Thus, 
below we discuss how the two factors, age of the stellar popula- 
tions and dust in the ISM and IGM, that contribute the most to 
deviations from a simple relation of proportionality between the 
A^f/y luminosity and the SFR affect these conclusions. 

The NUV is, in general, dominated by the contribution of 
stars with ages <100 Myr. The NUV images, which show patchy 
morphologies, give general support to this. But the median 
A'f/y profiles at z~0 are close to exponential, and show substan- 
tial differences from the B profile (a tracer of the distribution of 
significantly older stars, with ages > 1 Gyr), but not so "erratic" 
as typical profiles in Ha, a better tracer of on-going SF. In partic- 
ular, the central excesses of A^[/y emission that we find at z~0, 
as discussed above, are suspected to be, at least in part, "con- 
taminated" by the contribution of older stars (the signatures of 
bulges and pseudo-bulges). 

The dust is a more involved issue, that requires special atten- 
tion, as its absorption is more important at shorter wavelengths 
(see e.g. Boissier et al. 2004; Popescu et al. 120001 Tuffs et al. 
12004 ; Mollenhoff et al. |2006l l. We have used the prescription 
for dust attenuation as a function of radius in local disc galax- 
ies by Boissier et al. |2004| to correct the Bjif and NUVrf profiles 
for dust absorption, as is shown in Sec. 13.91 This dust absorp- 
tion is maximum in the innermost parts of the profiles, and is 
negligible beyond ~2.5-3 Reff , and so the underlying B^f and 
NUVrf profiles have presumably steeper inner profiles than ob- 
served. Furthermore, the absorption is higher in NUV ihan in 
B . We have explored how this dust-corrections affect the mea- 
sured parameters, and the stacked intensity and colour profiles. 
This test indicates several things. Part of the "truncation" effect 
in median A^ t/y profiles at ~/?(.// could be due to dust. The 
valley-like appearance of the colour profiles {NUV -B ) seems to 
be intrinsic to the distribution of stellar populations at z ~0, but 
at intermediate redshifts it dissipates, with the significant red- 
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dening in the central parts seemingly being due to the effect dust. 
Though the dust has a significant effect on most parameters, such 
as C42 , Reff , fraction of NUV flux within a fixed radius, it does 
not, on its own, seem to mimic evolution for some of the param- 
eters (e.g. Reff {NUV)), and for those which does, the effect is 
not great (e.g. C42 , fraction of flux within the central 2.5 kpc, 
Ayu), and largely compatible with the trends found before the 
corrections. Succinctly, it seems that the decrease in SFR sur- 
face density is more pronounced within the central 2.5 kpc, rel- 
ative to outer parts of the disc. This would explain the growth in 
Reff (NUVrf ) since z ~1 at a fixed mass, a phenomenon which 
is persistent after the dust corrections. Furthermore, and if we 
take these dust corrections for certain, it seems that the central 
excess of the Z? profiles, which we attribute to pseudo-bulges, 
could be largely in place at z ~ 1 , but it would be outshone by 
the more widely radially distributed emission from young stars 
at that time. This last interpretation neglects the effect of stellar 
mass redistribution within (dynamical evolution) and amongst 
galaxies (mergers) in the meantime, and so it could be erroneous. 

How can we understand the previous results in terms of the 
growth of stellar discs? We see that the (NUVrf -derived) SFR 
profiles are roughly exponential at z ~1, and since then the SF 
activity has decreased at all radii, but more so in the inner parts 
of the discs, leading to a downbending of the SFR profile with 
the knee at R ~3 kpc. In some sense, this is compatible with an 
inside-out scenario for the formation of the disc, but not exactly 
so. Up to a radius of R ~8 kpc it seems that there was more star- 
forming activity in the past, at all radii. But in the inner parts 
this activity has declined more rapidly since then. The net effect 
is that the stellar populations must have aged more in the central 
parts, and so they were formed before the outermost parts, which 
were being formed at z ~1, a process that is still going on nowa- 
days. In fact, we do not see that the SF has actually migrated 
to larger radii, within that radial range (imposed by the lower 
signal-to-noise ratio at large radii at z ~1). Instead, it seems that 
it is being been shut-off inside-out. Deeper images at all redshifts 
could provide valuable information with regard to what happens 
at larger radii. 

To end the discussion we comment on some observational 
biases, namely the change in spatial resolution with redshift, 
and the differences in the relative depths of the images. We have 
taken a careful approach to the resolution problem, and different 
tests (limiting distances to avoid the bias in the selection of the 
objects in the local sample, and comparing results obtained with 
images of equal, simulated, resolution) have shown that these 
results are robust in this regard. More delicate is the issue of 
the differences in depth, which are substantial, as commented 
in section [2.4.31 As we said there, the GALEX -NUV images 
are ~10 times deeper than the SDSS , and these are <20 times 
deeper than GOODS . So, the portrait here presented on the 
A^f/y profiles of disc galaxies in the past must be incomplete, 
as we do not reach the same depths as with the local data. But 
some of the results in Local and High-z sample are not consis- 
tent with absence of evolutionary effects (even when considering 
these differences in depth), such as the change in the profiles, 
with the development of substantial central emission at z~0, or 
the significant decrease in surface brightness at R^ff to a 20% of 
the value at z~ 1 . If there were significant, undetected, low sur- 
face brightness in the GOODS images this would make the in- 
tensity level at R^ff lower at z~l, and so the difference would be 
less significant, but that does not convey the fact that the z~l pro- 
files are brighter on average than those of galaxies at z~0, as the 
differences in the are not so big (R^y.y(z ~ 1) ~R^jj {z ~ 0) 
as to render these differences a mere selection effect. Still, the 



discs could have had larger R^ff 's in the past, as there could be 
a low surface brightness contribution that may be unnoticed at 
intermediate redshifts, but can be detected in local objects. For 
these reasons it would be desirable to have deeper rest-frame 
NUV and B images at z~l. 

5. Summary & Conclusions 

We have studied the radial distribution of flux in the rest-frame 
NUV and B bands for a sample of 270 disc galaxies, in the range 
0<z<l (i.e. exploring the last ~8 Gyr of cosmic age). The anal- 
ysis is focused on the evolution in the spatial extension (param- 
eterized by effective radii, Reff) and concentration (using ap- 
propriate parameters) of the rest-frame NUV flux, which may be 
taken as an indicator of recent star forming activity. In this way, 
and to our knowledge, this is the first systematic study of the 
radial distribution of star formation in disc galaxies extending 
through a cosmologically significant range of redshifts. Our aim 
is to provide further constraints on the mechanisms which shape 
the structure of stellar discs with epoch. These are our conclu- 
sions: 

1. At a fixed luminosity {M[,--21 mag) galaxies have larger 
values of Reff in the rest-frame NUV band, NUVrf , by a 
factor 1.5+0.1 at z =1 than at z =0 (a factor 1.4±0.1 in the 
rest-frame Z?band, Brf )■ 

2. At a given stellar mass (M*= 10'" Mq), the Reff of 
the galaxies in NUVRF^as slightly increased, by a factor 
1.18±0.06. The Reff in Bband, at a given stellar mass, has 
increased by a smaller factor, 1 . 10±0.05 in the same range of 
redshifts. This is in fair agreement with other observations, 
which report no evolution in that parameter at longer wave- 
lengths in this range of redshifts (e.g. Barden et al. 2005). 

3. We find no significant evolution in the concentration 
of A^f/y emission in this range of redshifts (using the 
C42 parameter, Kent |19851 l, nor in the fraction of flux con- 
tained in the central parts of the galaxies, when this central 
region is defined in relation to the Petrosian radius (in B ). 
Nonetheless, galaxies at lower redshift show a larger differ- 
ence in surface brightness, fj., from centre to Reff (i.e. an in- 
crease in this difference by ~0.8 mag in NUV and by ~1.6 
mag in Z? between z~0 and z~l). 

4. Galaxies at z~ 1 contain a larger fraction of NUV flux within 
R<2.5 kpc than those at z~0 (-30% vs. >15%). 

5. The averaged surface brightness profiles in rest-frame 
NUV and B bands for the galaxies under study show: a) the 
NUV and Z? profiles at z~l are quite similar, being somewhat 
"flattened" in the central parts, showing no evidence of a 
distinct bulge emission. At z~0 the profiles have a ~20% 
and ~40% of the brightness at z~ 1 in NUV and B , respec- 
tively (measured at the Brf effective radii in each case). At 
lower redshifts the profiles show clear signs of bulges or 
pseudo bulges, more prominently in the B band, which could 
be understood as a by-product of progressive accumulation 
of stellar mass in the inner parts of the galaxies. However, 
it could be that these bulges were partially outshone at z ~ 1 
by discs with their intense star forming activity. Analysis of 
dust-corrected NUV and B profiles seems to favour the sec- 
ond interpretation, but this is still inconclusive. 

6. At both redshifts, the median A^f/y profiles are similar to 
those of stellar discs which show truncations, with the cor- 
responding break roughly at the position of the (B -band) ef- 
fective radius. The median colour profiles (NUV -B) show 
an evident "blue valley" shape, with a minimum in colour 
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found around R^ff (Bsf ), a feature similar to that found in 
averaged colour profiles of truncated discs in galaxies at z ~ 1 
(Azzollini et al. 2008a) and in local galaxies (Bakos et al. 
I2008]) . 

7. The fraction of NUV fiux enclosed by the "break" radius, 
Rsr, in truncated galaxies is ~9Q% at z~0 and >60% (at 
0.5 < z <1.1), with significant dispersion. This finding seems 
to favour models for stellar disc truncations which postu- 
late that these are due to a cut-off in the star-forming ac- 
tivity at the "break" radius (e.g. Elmegreen & Parravano 
[19941 Scha ye 120041 Elmegreen & Hunter 2006 Debattista 
et al. 12006 ' Roskaret al. 2008a), though the lower fraction at 
higher redshift indicates that this interpretation may not be 
complete. 

Finally, it would be desirable to perform an analysis of stel- 
lar populations (e.g. from SED modelling) at different galacto- 
centric radii which would complement the study of the distribu- 
tion of NUV flux, in order to weight the effects of varying dust 
absorption, and the ageing of underlying stellar populations, and 
disentangle them from genuine evolutionary changes in the dis- 
tribution of star formation across stellar discs. 
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APPENDIX 

DEGRADATION OF IMAGES TO A COMMON SPATIAL 
RESOLUTION 

Here we explain the process by which we degrade the im- 
ages of galaxies in different redshift bins to a common 
resolution, in order to avoid potential bias due to resolution 
differences. These transformations are always degradations by 
definition, in the sense that we apply only transformations in 
which the final images have coarser sampling (larger angular 
scale) and a larger value of the FWHM than the initial images. 
Otherwise, we would be merely interpolating the images, which 
is not useful in this context, or would have to deconvolve the 
images from their PSFs, a process which usually generates 
artefacts in the images. 

First, we interpolate the pixel values of the GALEX- 
A^f/y images to change their angular scale y from 1.5"/pixel to 
the 0.396"/pixel value of the SDSS images. We checked that 
this does not produce significant alterations in the extracted ra- 
dial distributions of NUV flux. We then selected a projected di- 
mension of the FWHM in kpc, FWHMkpc , which is larger than 
the projected size of the FWHM of the GOODS -B435 images at 
z=l.l, and also larger than the projected size of the relatively 
large GALEX -NUV FWHM (5"), at the limiting distance of the 
Local sample, which is 60 Mpc. Once the FWHM^pc is fixed 
(=1.5 kpc), all images are convolved with a Gaussian kernel. 
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with a certain matching FWHM (F,,,), so as to produce an im- 
age whose final FWHM (Ff), subtends the desired FWHMkpc at 
the distance of the object. To compute this F^ we make use of 
the expression: 

F,n = ^jF]-Fl (.2) 

where F,, is the original value of the FWHM . This expres- 
sion is valid for Gaussian PSFs, which is our working assump- 
tion. Using expression and convolving the images accord- 
ing to the described scheme, we produce images of objects in the 
Local and High-z sample which share a common physical size 
of the PS¥,FWHMkpc- 

But the images still have different physical scales, F; i.e. their 
pixels still project different lengths, in kpc, at the distances of the 
targets. We select a certain value of F, that which corresponds to 
the 1.5" pixels of the GALEX images at a distance of 60 Mpc 
(r=0.5 kpc), which is also larger than the values of this param- 
eter for the galaxies in the High-z sample (F <0.25 kpc in that 
sample). Then, all the selected images are re-sampled to match 
that value of F by performing a suitable "binning" on them (i.e. 
averaging a certain number of pixels, though this number may 
be fractional, which is computationally feasible by means of in- 
terpolation). 

These procedures of convolution and re-sampling were per- 
formed using several IRAF tasks and a script coded in Python 
language. 

We have seen that there are also substantial differences with 
regard to photometric depth amongst the images from the sur- 
veys used. There is the possibility to adjust the intensities (scale 
them) and noise level (adding Gaussian nose) of the images of 
the galaxies in the Local sample to match the significance level 
of the images of the High-z sample. But this strategy has also its 
inconveniences, firstly a) the resulting images would not unam- 
biguously represent the appearance of local galaxies observed 
at High-z, given the uncertainties in the intensity scaling factor, 
through the corresponding k-correction (which is reasonable to 
assume as spatially variant); and b) even more important, the de- 
crease in surface brightness between z~l and z~0 has been so 
significant, that local objects would be, in many cases, practi- 
cally undetected in the GOODS-ACS images. This is why we 
have not attempted that type of simulation in this work. 
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